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1.0 INTRODUCTION 


In both the so-called Intensive and Global modes, the- GEOS-3 radar al- 
timeter acquires altitude and backscattered power data; furthermore, in the 
Intensive mode, instantaneous and/or average return waveform data are also^ 
obtained. Upon further processing, these data will yield measurements of the 
geo id undulation, the normalized surface scattering cross section (cr“), and 
the rms height of the small scale surface perturbations*. The key to the 
transformation from raw radar data to these- surface characteristics is an 
in depth understanding of the scattering of electromagnetic energy from 
rough surfaces. For most applications of the GEOS-3 data, the data trans- 
formation problem hinges, more specifically, on knowing what influences the 
time varying average return power waveform. 

The impact of the radar upon the post detection average return power 
waveform is contained in the three following system parameters; the video 
point target response, the tracking loop jitter density, and the pointing 
angle of the radar antenna electrical boresight with respect to nadir [1] . 

The point target response of the radar is measured by the Sample and Hold 
(S&H) gates during the RF BIAS TEST in the BIT/CAL mode of operation. 
Furthermore, the tracking loop jitter density can be determined by histo- 
graming the tracking loop jitter which is obtained directly from the track- 
ing loop output. Thus, the first two radar dependent factors are directly 
measurable and therefore can be properly accotinted for. The pointing angle, 
on the other hand, is not measured and only its estimated range of variation 
is known. For certain early-mission orbits, estimates of the spacecraft 
z-axis pointing with respect to nadir were computed [2] . However, since the 
angle between the spacecraft z-axis and the radar antenna electrical bore- 
sight is not precisely known, such results do not help matters much. 

Since an uncertainty in the actual pointing angle of the radar antenna 
can lead to larger error bounds on radar derived estimates of surface char- 
acteristics, it is clearly in order to investigate means for determining the 
pointing angle. This report presents the results of a study to estimate the 
radar pointing angle through the use of the average Plateau gate output and 


*Small scale surface perturbations, as used here, refer to those surface 
features whose spatial wavelength is much smaller than the effective spot 
size of the radar. 



the average Attitude/Speciilar gate output. In addition, error bounds on the 
estimate are obtained based upon the number of returns used to form the aver- 
age. The pointing angle induced altitude bias is computed for the Global 
mode* end the variation of a“ with pointing angle is also presented. 

1.1 Summary of Results 

This report derives a relationship between the ratio of the outputs of 
the Attitude/ Specular and Plateau Integrating gates and the angle of the 
radar antenna boresight relative to nadir. The relationship is shown to be 
single valued and therefore amenable to use as a means of estimating the 
radar pointing angle. Estimation curves are obtained for both the Intensive 
and Global inodes. Stationary approximations to the autocorrelation function 
for the predetection backscattered signal are obtained which are applicable 
' over the delay extent of the Plateau and Attitude/ Specular gates. With 
the inclusion of receiver noise and rather simple models of the IF and video 
filters, estimates of the ten second standard deviation of the pointing angle 
estimate are obtained as a function of the effective receiver signal-to- 
noise ratio. The one-sigma statistical error of the method is on the order 
of 0.33° at nadir for Global Mode and 0.23° at nadir for the Intensive Mode. 
This error decreases to less than 0.07° for both modes when the pointing 
angle exceeds 1° . Based upon the statistical error of the technique, the 
method is best suited for use when the pointing error is on the order of 
one half a beamwidth. 

The effect of inadequate temperature dependent calibration of the 
Attitude/Specular and Plateau gates is shown to be a nonnegligible source of 
error in the estimation process. A five percent change in gate calibration 
is shown to give rise to an error which is a nonlinear function of the true 
pointing angle, i.e. the error is largest at nadir (»0.4°) and decreases to 
about 0.1° at one half a beamwidth pointing error. This error source is 
hypothesized to be the cause of disagreement between preliminary pointing 
angle estimates obtained during Global Mode and spacecraft-based sensor es- 
timates . 


*Attitude induced biases have previously been presented for the Intensive 
Mode in [1]. 
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Curves are obtained for converting received ' power measurements into 
pointing angle corrected estimates of <j® for both Global and Intensive Modes. 
Waveform and attitude induced altitude bias is confuted for the Global Mode 
and shown to be a major source of error relative to absolute altimetry. The 
effect of pointing angle upon waveheight estimation is discussed and dis- 
missed assuming that the pointing error is less than one half a beamwidth. 

The problem of a discrepancy between the AGC Clean and Clutter Inten- 
sive Mode calibration data is addressed and resolved using an analyti- 
cal approach. The sample-by-sample AGC control voltage is modeled as an 
average of logarithms of the Plateau gate outputs. Because of this loga- 
rithmic dependence, the average of the AGC control voltage differs from 
the deterministic (Glean) result by a constant value, i.e. a bias. This 
bias is equivalent to .a 2.5 dB separation between the Clean and Clutter 
calibration curves. The actual data indicates a separation of 3.5 to 4 dB; 
thus, there is a 1 to 1.5 dB residual error which cannot be accounted for. 
Considering the complexity of simulating a noise-like coded return, this 
residual error may well be due to experimental error. To within this re- 
sidual error, the Clutter AGC curves are shown to be the proper calibra- 
tion data for reducing inflight data. The bias for Global Mode is shown 
to be negligible, because the bias depends inversely upon the ntanber of 
returns averaged prior to entering the AGC loop. For -the Global Mode, the 
number of returns averaged was -sixteen and the resulting bias is shown to 
be -0.14 dB. The analysis presented here not only resolves the GEOS-3 
Clean vs. Clutter problem, but also sheds more light on the operation of 
logarithmic AGC systems for extended target scattering. 

The final two sections of the report deal with the implementation and 
-evalxiation of the models developed in the previous sections. A simple 
closed-form algorithm for converting the angle estimation function into 
pointing angle estimates is presented for both operating modes. A limited 
comparison of radar derived angle estimates with attitude control system 
results indicates the apparent presence of a 0.6® to 0.8° bias in the Global 
Mode. This bias is attributed to inadequate integrating gate calibration. 
Intensive Model comparisons indicate reasonable agreement with no apparent 
bias. A tjrpographlcal error in the Plateau gate calibration data has, 
unfortunately, been Implemented in the basic data correction/calibration ' 
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program at WFC. Temperature dependent correction routines have been de- 
veloped to overcome this problem in an hoc fashion. 

Closed-form algorithms have also been developed for computing from 
return power measurements. These algorithms also depend upon transmitted 
power and pointing angle. Global Mode cr“ values compare very favorably 
with Sky lab results for similar surface conditions. Intensive Mode values 
are some 3 dB lower than the Global Mode results. This was attributed to 
the use of Clean rather than Clutter AGC calibration data. An in-orbit 
test was conducted and it verified that the Clutter calibration data is the 
proper set. Correction routines have been developed to overcome the errone- 
ous data calibration problem. 

2.0 ESTIMATION OF THE RADAR POINTING ANGLE 

2.1 Background 

During study of the Skylab radar altimeter data, it was found that the 
pointing angle of the radar antenna boresight relative to nadir could be ac- 
curately determined from the shape of the post detection average return 
waveform [3].' This finding was a direct consequence of the fact that for 
near nadir pointing only the trailing edge portion of the AGC normalized 
average return was sensitive to a change in pointing angle. Figure 1 illus- 
trates the behavior of the Skylab average return waveform for two different 
pointing angles. Of particular note in this figure is the fact that the 
leading edge of the average return is relatively insensitive to the change 
in pointing angle while the trailing edge exhibits a marked change. The 
pointing angle estimation procedures developed for Skylab involved compari- 
son of the measured post detection average return waveform (in the trailing 
edge portion of the return) with computed waveforms- [7]. Such a direct com- 
parison was possible for Skylab because waveform samples of the return were 
obtained well into the trailing edge portion of the return. For a 1500 
pulse average, the one-sigma error in the pointing angle estimate due to 
the statistical nature of the return was shown to be less than 0.05® 17] 
for, roughly, one half a beamwldth pointing error. 

For the GEOS radar altimeter, no waveform point samplers, i.e.. Sample 
and Hold gates, were located sufficiently far into the trailing edge portion 
of the return to directly apply the estimation procedures developed for Skylab. 
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However, the General Electric Company originally proposed an alternate im- 
plejnentation [4] whereby an Integrating gate (called the Attitude/ Specular 
gate) was located in the trailing edge portion of the return; the location 
of this integrating gate relative to the Ramp and Plateau (tracking) gates 
is shown in Figure 2 for both the Intensive and Global Modes, By comparing 
the output of the Attitude/Specular gate (which is' sensitive to the pointing 
angle) to the output of the Plateau gate (which is much less sensitive to 
the pointing angle) the radar pointing angle could be Indirectly measured. 
Preliminary analyses of this approach [4,1,5] indicated that the technique 
showed sufficient promise as to warrant further in depth analysis. The 
purpose of this section is to present the analysis for determining the es- 
timated pointing angle and its statistical precision from the Plateau and 
Attitude/ Specular integrating gate- outputs. 

2.2 The Estimation Function 

Figure 3 is a block diagram [6] of those parts of the GEOS radar altim- 
eter receiver which are essential to the ai^le estimation problem. It is 
important to note that the outputs of the integrating gates are passes 
through a one second averager before being telemetered to earth. The first 
step in the analysis comprises relating the filtered outputs of the inte- 
grating gates to the radar pointing angle. 

2.2.1 The Average Integrating Gate Outputs 

For the ith return voltage, the output of the IF filter/amplifier, 
which also contains the piilse compression network for the Intensive Mode 
(IM), may be represented as 

x^(t) * l^v^^(t) +n^^(t)] cosM^t - sinw^t , (1) 

where is the IF center frequency. The in-phase and quadrature return 
voltage components v . (t) and v . (t) are independent, zero mean, Gaussian 

Cj SJ ^ 

random variables with a time varying variance equal to KF^(t), i.e. the 

product of the average return power and a constant which depends upon the 

RF gain of the receiver and the A6C. Similarly, the in-phase and quadrature 

noise voltage components n , (t) and n ,(t) are Independent, zero mean, 

c^ S’! 
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TRACKING GATES 


INTENSIVE MODE 

Tq = 1 2.5 ns. 

T| =62.5 ns. 

T£ =700 ns. 
1 ^ 3=200 ns. 



RAMP 


PLATEAU 



A, 

ATTITUDE/ 1 
SPEC. 1 

Tr 



Y 



T. 



^ T 




>08 ^ 

1 ^ i 



GLOBAL MODE 

T (3 = 200 ns. 
Tj'=300 ns. 

Tg^TOp ns. 
Tgg=200 ns. 


Figure 2. Nominal GEOS tracking gates' lengths and locations 
for both Intensive and Global Modes. 



Figure 3. 


A simplified block diagram of the GEOS-3 radar altimeter 
receiver pertinent to the attitude estimation process - 
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Gaussian random variables xfith variance KN where N = kT B„_ and k is 

o » o e IF 

Boltzman's constant, is the effective receiver noise temperature and 
Bjp is the IF filter/amplifier noise bandwidth* Squaring (1) and regrouping 
terms yields, for the output of the square law detector, 

='3 (« - i I [vej (« (t)] ^ + [»s3 +”s3 («] I 

2 ■^“c3 ■ ['s3 '^“s3 ^ I 

- <« +“<.j c«j [vgj (« +”sj («] sln2u^t (2) 

Since the video filter/amplifier has a low-pass characteristic with a band- 
width, much less than 2f , its output is given by 

( 

^3 ■ 2 I [''c3 ■^“c3 <*>] ^ ^3 ■^“s3 ^ I “v'*’ ’ ■ 


where h^(t) is the impulse response of the video filter /amplifier and the 
symbol fi denotes convolution. The series capacitor between the video ampli- 
fier and the two gates essentially blocks those components of (t) which, 
in the mean, are either constant or slowly varying with time. Thus, the 
input to the two gates is given by the following; 

= Zj(t) - (t)] (« 

2 

since n (t)0h^(t) is the only term which yields a constant mean value. 

The average output of either the Plateau or Attitude/ Specular integrat- 
ing gates is given by 


T + 6 


e •= E{e(j)} = G ^ E|z^^(t) | hg(T + 6 -t)dt 


(5) 


where T is the time (in two-way time delay coordinates) at which the particu- 
lar gate is switched on, 6 is the gate length, G is the gain of the gate, and 
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hg(t) is the norsaalized impulse response of the integrating gate. Using (3) 
and (4), the average return voltage into the integrating gates is related to 
the IF average return power as follows; 

E|z^^(t)| = KP^Ct) ©h^(t) (6) 

The IF average return power may be expressed [1] as a convolution of the 
IF point target power .response*, P^pCt), with the so-called flat surface 
impulse response, P^gCt). Interchanging the order of convolutions, (6) may 
be rewritten as 

where V^(t) is the convolution of with the video filter /amplifier im- 

pulse response. V^(t) is the video voltage point target response and it is 
equivalent to _the voltage waveform at the output of the video filter/amplifier 
during the so-called RF Bias Test of BIT/CAL for the GEOS altimeter. Sub- 
stituting (7) into (5) yields the following expression for the average out- 
put of either the Plateau or Attitude/Specular gates; 


T+6 


e = E|e(j)} ^ V^(-t)«Ppg(t)hg(T+6-t)dt 


( 8 ) 


According to [6], the integrating gates are analog integrating circuits 
with a time constant equal to four times the gate width; thus, the normalized 
impulse response of the integrating gate is given by 


hg(T+5-t) 


exp [^(T + d - 1) /4d] 
yexp [(T + 6 - 1) /4d] dt 


(9) 


*For the Intensive Mode, the IF point target response is equivalent to the 
squared ambiguity function of the transmitted signal. For the Global Mode, 
the equivalence is with the cross-ambiguity function squared, because matched 
filtering is not used. 
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For near nadir operation, the appropriate flat surface impulse response 
has been shown [1] to be given by the following; 




G X c )//<>/ \ 

o o fh.2f.Uc. 

9 ? exp ( - - sin § - “ t cos2? ) 

4(4tt)^L L h-^ \ Y ^ / 

^ ' p r 


( 10 ) 


for t > 0 and ^j.g(t) = 0 for t < 0. In (10), is the peak transmitted 
power, Tij.p is the measured pulse compression ratio, 6^ is the boresight 
gain of the radar antenna, X is the radar wavelength, c is the speed of 
light (0.3m/ns.), o® is the surface scattering cross section per xmit scat- 
tering area, is the two-way path attenuation, comprises the appro- 
priate radar losses, h is the altitude of the radar above the mean surface, 
and Y is a parameter related to the beamwidth of the antenna pattern (1.49><10 
for GEOS-3). Also, in (10) and the remainder of the material in this report, 
time is measured relative to the total two-way delay time (2h/c) from the 
radar to the mean surface. The variable 5 in (10) is the pointing angle; it 
is defined as the’ angle between the electrical boresight of the radar an- 
tenna and the nadir axis. The angle is defined by the following relation- 
ship; 




1/2 


( 11 ) 


The video voltage point target was measured during preflight testing 
of the altimeter [6-], and the results are shown in Figures 4 and 5 for the 
Intensive and Global Modes, respectively. A discretely sampled version of 
Figure 4 is also available from BIT/CAL data (KF. Bias Test) obtained during 
actual operation of the altimeter. Examination of RF Bias Test data [8] ob- 
tained during Inflight operation of the altimeter has indicated no signifi- 
cant changes in the Intensive Mode video point target response. No record- 
ing of the video point target response is accomplished in BIT/CAL for the 
Global Mode. 
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NORMALIZED AMPLITUDE 


LO - 
0 . 8 - 
0 , 6 - 

0 . 4 - 
0 . 2 - 
0 - 



2.2.2 The Qae Second Hardware Averager 


The output of either the Plateau or Attitude/Specular integrating gate 
is, according to [6] and Figure 3, averaged for one second and then sampled 
by the telemetry system. The averaging process is accomplished by an ana- 
log filter, however, the transfer characteristics of this filter were not 
specifically detailed in [6]. For this reason, the analog averager or fil- 
ter is modeled by a discrete process in which the output Is formed by aver- 
aging the number of returns In a one second time Interval. That is, the 
output of the one second averager in Figure 3 is approximated by the follow- 



( 12 ) 


where M is the number of returns integrated by either the Plateau of Atti- 
tude/Specular gates in a one second interval. For the Intensive Mode, M 
is equal to 100. In the Global Mode, M = 1600 for the Plateau gate and 
M ® 100 for the Attitude/Spectdar gate; the difference is dtie to the fact 
that only the first rettim of the sixteen pulse hurst (per prf ) ' is inte- 
grated by the Attitude/Specular gate [6]. Table I summarizes the important 


TABIE I 


A Tabulation of Intensive and Global Mode Plateau 
and Attitude/Specular Gate Parameters 


GLOBAL MODE 


INTENSIVE MODE 


Plateau Attitude/ 
Gate Specular Gate 


'as(g) 


Plateau 
' Gate 




Attitude/ 
Specular Gate 


*'as(g) 


300 ns 


200 ns 


700 ns 


200 ns 


62.5 ns 


700 ns 


Relative Gate 
Tum-On Time, T 


Gate Width, 6 


Number of Pulses 
Avgd. lu 1 Sec, M 


*These gains are functions of the tracker temperature and input signal 
strength (near saturation); extensive calibration curves are given in 
[6]. The same Attitude/ Specular integrating gate is used for both 
Global and Intensive Modes. 
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gate parameters for the Plateau and Attitude/ Specular integrating gates for 
both the Global and Intensive Modes. 

2.2.3 The Minimum Averaging Time and Its Effect Upon the Estimation Function 

The ensemble average of the one second averaged integrating gate out- 
put is equal to the average gate output, >1.6., e = e(i). Using equations 
(8), (9), (10) and the waveforms in Figures 4 and 5 along with the gate pa- 
rameters in Table I, the average output of the gates can be computed. ‘ It 
is now possible to construct a function of the Plateau and Attitude/Specular 
gate outputs which will result in a single-valued determination of the pointing 
angle; this function will be called the estimation function. It is obvious 
from (8) that this estimation function should be dependent upon the ratio 
of the Plateau and Attitude/Specular gate outputs in order to eliminate the 
unknown AGO dependent constant K. Furthermore, the estimation function 
•should also be a function of the ratio of the gate outputs to the gate gains 
since this eliminates the need for detailed temperature corrections and is 
also consistent with Wallops Flight Center data processing products. A 
previous study has suggested the following form [1] for the estimation func- 
tion; 


A = 1 - 


e /G 

as as 

e /G 

, P P 


(13) 


where and are the one second averaged outputs of the Attitude/Specu- 
lar and Plateau gates, respectively, and G and G are the gate gains. 
However, defining the estimation function in terms of the one second aver- 
aged outputs, as above, suffers the drawback that the average estimation 
function. A, does not converge to the ratio of the average gate outputs. 

That is, since c and e are derived from at least a 100-pulse average, 
p as 

they will be essentially Gaussian with the following means and variances; 


e = e 
as as 


e = e 

P P 


( 14 ) 




Var(e ) 


as 


M 


Var(e ) 

Var(e ) :s ^ 

P 


as 


M 


(15) 


IS 



where M and M are the number of pulses averaged to form the one second 
as p 

averaged Attitude/ Specular and Plateau gate outputs (see Table I). In 
order to determine the average estimation function, the joint density func- 
tion of e G /e G must be computed. This density function can be obtain- 
as p p as 

ed by the methods given in [9]; however, because the two processes e and 

P 

e are nonzero mean Gaussian, the first moment is not easily calculated, 
as 

An alternate approach [9, pg. 212] to computing the average of the estima- 
tion function is to expand (1-e /e G ) in a Taylor series about 

as p p as 

e ~ and e = T and only retain the significant terms (after aver- 
as as P P 

aging) . This procedure is valid only when the probability masses of 

and e are very concentrated near their "center of gravity", i.e.,e and 
p as 

, and (1 - s G /e G ) is smooth in the vicinity of this point. Refer- 
p as p p as 

ence 10 indicates the order of error this approximation can lead to when 
the above assumptions are violated; however, because of the degree of vari- 
ance reduction brought about by the one second averaging, this approach can 
be safely used to compute the average estimation function. Therefore, using 
the Taylor series approach, the average estimation fxinction can be shown to 
reduce to the following form; 


e G 

1 

e G 
p as 


Var(e^) 
1 + 2 - 


/- x2 


Rewriting (16) in terms of the average gate outputs yields 


A = 1 - 


e G 
as p 

e G 
p as 


1 + 


Var(e ) 
P 

M (e )^ 
P P 


(16) 


(17) 


Equation (17) illustrates why it is not desirable to define the esti- 
mation as in (13); that is, the average estimation function depends upon 
the variance of the Plateau gate output. For the Global Mode, this depen- 
dence does not create a problem since M is 1600 and Var(e ) is relatively 

p P 

Insensitive to receiver noise. However, for the Intensive Mode, M = 100 

P 

and Var(e ) is more strongly dependent upon the receiver noise because of 
P 

the lower signal-to-noise margin- Thus, if the estimation function is de- 
fined as in (13) , the average estimation function will depend upon the 
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receiver signal-to-nolse ratio and this is not a desirable situation. The 
obvious way to overcome this problem is to increase in (17) by at least 
an order of magnitude; then, the variance dependent term in (17) will be 
negligibly small aiKi can be ignored. For this reason, the estimation func- 
tion A should be defined in terms of at least a ten second average of the 
gate outputs, i.e . 


where 


e /G 

= X - 

e /G 
P P 


lOM 


as 


e 

as 


ICM 


■z 

j=l 


e (j) 

as 


e 

P 



(18) 


(19) 


( 20 ) 


The remainder of this report will deal exclusively with the estimation fimc- 
tion as defined in (18). 

Ignoring the variance dependent term in (17 ) , the average estimation 
function becomes 


A = 


1 - 


e 

as 



e /G 
P P 


( 21 ) 


In (21) , the gain normalized gate outputs are directly compatible with the 
gate output data products provided by Wallops Flight Center. Figure 6 illus- 
trates how the average estimation function depends upon the radar antenna 
pointing antenna for both the Global and Intensive Modes. The greater slope 
of the curves as the pointing angle approaches a beamwidth (2.6®) reflects 
the fact that the average return waveform in the vicinity of the Attitude/- 
Specular gate changes more rapidly as the pointing angle approaches a beam- 

width. The smaller slope of the Global Mode curve is a direct consequence 
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AVERAGE ESTIMATION FUNCTION 



Figure 6* The average estimation function for Global and Intensive Modes* 


of the placement of the Plateau gate on the average return by the tracking 
loop time discriminator. That is, for the Global Mode, the center of the 
Plateau Gate occurs approximately 400 ns after the start of the return (see 
Figure 2)} at this relative delay time, the return shows a significant 
"droop” due to antenna pattern effects. Thus, both the Plateau and Attitude/- 
Specular gates encounter portions of the return which are dominated by an- 
tenna pattern induced droop. Similarly, as the pointing angle increases 
over the range shown in Figure 6, both gate outputs will also increase but 
at a different rates. For the Intensive Mode, the Plateau gate is located 
much nearer zero relative delay time and therefore is not effected by the 
antenna pattern; thus, the curve in Figure 6 for the Intensive Mode is pri- 
marily determined by the waveform in the vicinity of the Attitude/ Specular 
gate . 

2.3 The Precision of the Estimate 

. The previous section of this report was concerned with obtaining a 
relationship between the radar antenna pointing angle and the average Pla- 
teau and Attitude/ Specular integrating gate outputs. This section will be 
devoted to computing the statistical uncertainty of the pointing angle esti- 
mate. This uncertainty is a direct consequence of the fading and fluctua- 
ting nature of the i;etura and the nonnegligible influence of receiver noise. 
Since the analysis will be concerned with the variance of the Integrating 
gate outputs, this phase of the investigation is closely related to the 
problem of determining the precision of the GEOS-3 split-gate tracker de- 
rived altitude estimates [11]. 

2.3.1 The Variance of the Gate Outputs 

The average output of the integrating gates is given by (5); the vari- 
ance of the output is as follows; 

T+5 T+6 

Varte) = Gy f 
T T 

“ ( 22 ) 

IQ 



Siibstitutiag (3) and (4) In (22) and using the independent Gaussian prop- 
erty of v(t) and n(t) to simplify the resulting expectation operations [12], 
(22) reduces to the following form* ** 


T+6 T+6 


Var (e) 


- y*| |^R^(t^,t2) + 2R^(t^,t2)R^(t^,t2> + 


©h^(t^)®h^(t2) j hg(T+6-tj^)hg(T+5-t2)dt^dt^ (23) 

where R (t, ,t_) and R (t, ,t„) are the autocorrelation functions for the IF 
viz n 1 z 

in-phase or quadrature components of the return voltage and the noise, re- 
spectively. That is, the iatrapulse autocorrelation function for the re- 
turn is 

while the noise autocorrelation function is 


2.3.2 The Return Signal Autocorrelation Fxmction 

The autocorrelation of the noise is completely specified by the equi- 
valent receiver noise ten^erature and the bandpass characteristics of the 
IF filter /amplifier. The return voltage intrapulse autocorrelation func- 
tion is considerably more complex because, in general, it is nonstationary. 
Berger [13] has obtained an analytical expression for R^(t^,t 2 ) under cer- 
tain restrictive conditions. Generalizing his result and correcting some 
of his algebraic errors, it can be shown that the IF return voltage auto- 
correlation function is given by the following , 


*Itt the transformation from (22) to (23), it was assumed that the video 
filter /amp did not appreciably alter the shape of P^(t); measurements pre- 
sented in [6] indicate that this is a reasonable asstimption. 

**Equation (24) is valid if the scattering surface is noncoherent and the 
doppler spread in the return is small relative to the doppler extent of the 
radar ambiguity function; both of .these conditions are satisfied for over- 
ocean operation of the GEOS-3 radar altimeter. 
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( 24 ) 


f 1/2 1/2 

R^(tj^»t2) a; K / (t^- t) -t) Ppg(t) dt 

O 

1/2 

where P^^ (•) is the square root of the IF point target power response and 

Ppg(.) is the flat surface impulse response^ see equation (10). Equations 
(24) and (23) indicate a five-fold integration to accomplish the transfor- 
mation from known quantities, i.e., Pj^(t) and Pyg(t), to the variance of 
the gate outputs. For this reason it is most necessary to accomplish as 
many of these integrations in closed form as possible. 

2. 3. 2.1 Intensive Mode 

Inspection of the IF point target power response measurements present- 
ed in [6] shows that a Gaussian fom is a very good approximation for 
1/2 

Pjj, (t^) in the Intensive Mode, i. e. 

where a. z 5.3 ns and t z 2}/2 a. is time shift introduced to insure the 
i o i 

proper timing between P^ and Pj,g in (24). In (24), since Ppg(t) is nearly 

constant over the range of integration for which the product 
1/2 1/2 

^IF ^^l”*'^^IF ^^ 2 ”*^^ nonzero, the following approximation is valid; 


^^(t^.ta) 2: 






Substituting (25) in (26) yields 

R^(ti,t2) = • "^o) ] 


f 

1 

1 

"Vt - 2t 1 

i + 

— erf 

X ^ 0 

2 

2 

- 2^/z -* 


(26) 


(27) 


where erf(*) is the error function. Those parts of R^(t^,t 2 ) which give 
rise to its nonstationary nature are functionally dependent upon the sum 
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of the two time coordinates t^ and t^. As noted previously, the flat sur- 
face impulse response varies rather slowly as t^ and t 2 range over the ap- 
propiate time extent of the integrating gate, i.e. from T to T + 6. Thus, 
the flat surface impulse response may be adequately approximated by its 
value at the mid-point of the gate, i.e. 


FS 





- cj 


(28) 


Using the values for T and 6 shown in Table I, it may be verified that the 
argument of the error function is greater than two for both the Plateau 
and Attitude/ Specular gates; thus, since erf(t^2) ~ 1, 


erf 


V^2-"^o 

2i/2 a. 


1 


9 


and (27) simplifies to the following; 


R^(t^,t2) « K / 27ra^^ Ppg(l + 6/2 - t^) exp|^-(t^-t2)^/8a^J (29) 

Thus, for the Intensive Mode, the predetection intrapulse autocorrelation 
function is essentially stationary over the time extent of the Plateau and 
Attitude/ Spectilar gates. 


2. 3.2. 2 Global Mode 

For the Global Mode, the situation is more complex because the shape 
of the IF point target power response is not easily represented by a simple 
function (see Figure 77 of [6]) such as a single Gaussian as in the case of 
the Intensive Mode. To take advantage of some of the approximations de- 
veloped for the Intensive Mode, the square root of IF point target response 
was approximated by a four term series of Gaussian functions, i.e. 

4 

^I^^^t) - S exp [- ] (30) 

The amplitudes (a ) , time shifts (3 ) , and spread factors (O ) were generated 
mm m 
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by a computer program which minimizes the mean squared error between the 
measured function and the approximating series. A tabulation of these 
factors is shown in Table II. 

Substituting (30) in (24), xising the same rationale to remove the flat 
surface impulse response from inside the integral as in (26) , and integrat- 
ing term by term yields the following result; 




k=l m=l 


Fs( 


<V^m> 


cy^(t2- 


2 2 

k m 


-)• 


a a; 0 0- 
“ K m k 


m 




+ a. 


TABLE II 


Amplitudes, Time Shifts, And Spread Factors For The 

Four Term Gaussian Approxi m ation To The Square 

Root Of The IF Power Point Target Response (Global Mode) . 


m 

a 

m 

m 

0 

m 

1 

0.615 

55 

25 

2 

0.705 

125 

50 

3 

0.412 

190 

50 

4 

0.66 

300 

85 


Note: the units of 3 and a are nanoseconds 

m m 
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• exp 


(^1 ^ 2 ^ 


1 

2 


+ ~ erf 


( V^.> ( V^k>' 


/2 a, 


r~2 T 

a va + a 


k m m 


( 31 ) 


Using the parameter values in Table II , it can be shown that the argument 
of the error function in (31) is always greater than two as t^ and range 
over the Attitude/ Specular gate. When t^ and t^ range over the Plateau 
gate, the argument is also greater than two except for pairs of indices 


k = 3 , m = 4 

k=4 , m=3 

and k = 4 , m = 4 

Even for these terms, the error ftmction argument is not appreciably small- 
er than two; also, since there are only three out of the sixteen terms in 
(31) for which the error function argument is less than two, their effect 
upon the sum will be further diminished. Consequently, the error ftmction 
in (31) will be taken as one for both the Plateau and Attitude/Specular 
gates. 

A remaining problem with (31) is that the argument of the flat surface 
impulse response depends upon both summation indices and it cannot be re- 
moved from inside the double sum. However, because of the exponential fac- 
tor in (31) and the slowly varying nature of the flat surface impulse re- 
sponse, each term in the double series will be significant only in the 

neighborhood of t, - t- = 3, - B or t_ = t„ + 3, ~ 3 • With this transforma- 
liKm 12km 

tion, the argument of Ppg becomes t^ - 3j^ and (31) reduces to 


R (t,,t.) 

V 1 2 


k=l m=l 


a a, a a. 
m Tc m k 


^ / g ^+ O, 


m 


2 

k 


• exp 






(32) 
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The in the argvment of P^g 


is now replaced by an average value given by 





= 167.5 ns 


and removed from Inside the summation. This step is once 

again justified by noting that the flat surface ia^ulse response is rela- 
tively slowly var 3 ^g with time. As in the case of the Intensive Mode, t« 

^ 2 

in the argument of P^g is replaced by T +|- , i.e. the location of the mid 
point of either .the Plateau or Attitude/ Specular gate. Equation (32) there- 
fore reduces to the following form; 

+ 


= K.®r Pj/i 



• e^ 




(33) 


It should be noted that (33) is a symmetric function of t^ - because of 
the double summation. 


The double sum in (33) was computed as a fimctlon of the time differ- 
ence (tj^ - 1^) and compared to a single Gaussian fmction. The following 
approximation was found to be very good; 


tt 


a a, a a, 

m K m k 



exp 


[(V^2)-(gA)]" | 

+ "J) j 


(34) 


where O -102 ns and a - =161.37 ns. Figure 7 shows the comparison be- 
go gl 

tween the normalized double sum and the single Gaussian approximation. The 
final form for the Global Mode autocorrelation function is 
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NORMALIZED AMPLITUDE 



Figure 7 . Comparison of signal autocorrelation coefficient 
and a Gaussian approximation for Global Wtode. 
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( 35 ) 


where 

a 
go 

3 

1 

gl 

It should be noted that the end purpose of the above sequence of approxima- 
tions was to obtain an integrable form for the autocorrelation func- 
tion. In performing these approximations, the slight nonstationary nature 
of the autocorrelation function in the Plateau gate region has been Ignored 
and the resultant approximation is, no* doubt, a smoothed version of the 
true function. However, equation (35) does contain the basic form of the 
Global Mode autocorrelation function and is probably more than adequate for 
a variance calculation. 

2.3.3 The Noise Autocorrelation Ftmction and Receiver Model 

Now that the signal autocorrelation functions are known, the noise 
autocorrelation fimction must be obtained. As [11] Indicates, even for a 
simple model of the IF and video filters, the problem is extremely tedious. 
To simplify the labor as much as possible and yet retain the basic character 
of the problem, the IF filter was modeled as a two-pole RC filter while the 
video filter was taken to be a single-pole RC filter. The 3 dB bandwidths 
were matched to the measurements; presented in [6]; namely, for the IF filter 

Intensive Mode: ~ 52.4 MHz 

Global Mode : BW^^, = 40.4 MHz 

and for the video filter 

Intensive Mode: BW = 50 MHz’ 

v 

Global Mode : BW^ = 5 MHz 


= 102 ns 

» 167.5 ns 
= 161.37 ns 
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For the two-pole IF filter, the autocorrelation function of the noise is 
given by the following; 


R^(ti,t2) = K(1.22kT^)(BWjp) |(4.88BW^)|tj^-t2| +l| 
• exp 4 , 88 BWjp j t^“t 2 1 ] 


(36) 


where k is Boltzman's constant and is the equivalent noise temperature 
of the receiver front end. » 


Before proceeding further with the determination of the variance of 
the output of the integrating gates, equation (23) will be further simpli- 
fied. First, in the double convolution of the video impulse response with 
the product of the noise and signal autocorrelations, the following approxi- 
mation is made; 


R,(ti.t2)E„(ti,t2)9h^(ti)9h^(t2) = R^(t2.t2)[E^(t2,t2)9h^(t2)9h^(t2)] (37: 


Equation (37) assvimes that the noise decorrelates much more rapidly than the 
signal and, therefore, the video filter essentially only effects the noise. 
For the Global Mode, this is the case as will be shown later in this section. 
For the Intensive Mode, the video filter has little effect upon either the 
signal or the noise; thus (37) is more in the form of an identity, as will 
be shown. Also in equation (23) , the exponential weighting of the integrat- 
ing gates will be simplified to a uniform weighting, i.e. for t In (T,T + cS) 

hg(T + 6-t)::| (38) 


This simplification may be justified by noting that since the autocorrela- 
tion functions in (23) decorrelate rather rapidly over the range of inte- 
gration, i.e. the extent of the gates, the integrations in (23) are insen- 
sitive to minor changes in the gate weighting. In other words, the inte- 
grals in (23) are much more sensitive to the autocorrelation functions 
than the form of. the integrating gate impulse response. Thus, using the 
above approximations and the fact that the autocorrelation functions are 
essentially stationary, equation (23) may be simplified to the following 
form; 
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u 


where x has, replaced and 


R^(T) = R^(T)0h^(T)0h^(x) 


R^(T) 

n 


R^(x)®h (x)6h (X) 


R,(T) 


\(T)0h^(T)®h^(T) 


For the noise terms in (39), the double convolutions may be accomplish- 
ed either directly or in the spectral domain. The calculations are straight- 
forward but extremely tedious, as noted previously, so only the results will 

( 

be presented. For the Global Mode, 


-aalTl 


R^(x) = K^fo.SlkT 1.1845 e -0.171 e 

n L e IFJ 


-2a^ T 


-2ajx| 


- 1.518(a2|x|) e 


+ 1.455(a,lx|) e 


-2a^ T 


-a- jx 


Rj^(x) »= 0.618 k[o' 51 RTeBj-p] 1 1.319 e - 0.319 e 


-2a^ x| 


- 0.642(a„|xl) e 


- 2aj^|xj , 


where a. = 2 ttB and a, « 4,88 B,^„. For the Intensive Mode, 
2 V 1 IF 
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-021x1 - 2a2^1x 


R^(t) = j 8.535 e -7.535 e 


-2aJx- 


2 . 


2a. (t 


- 2.532(a^jxj) e - 0.395(a^|x|)'^ e 


( 42 ) 


and 


\ (x) = 0. 964 K [l. 01 k B J 


-a^K 


-2aj^lxj 


11.296 e 


- 10.296 e 


+ 2.897(0 [x|) e 


- 203 ^ M 


(43) 


The effective receiver noise power at the receiver Input, referred to the 
effective system bandwidth, is equal to /K • 

2.3.4 Comparison of Autocorrelation Coefficients 

From eqtjatlons (29) and (35), the square of the return autocorrela- 
tion function may be written in the following form; 


R^(x) 

V 


= R^(0) exp [•- x^/2a^] 


(44) 


where, for the Global Mode, 


r2<0) 




0=0 


gl 


(45) 


and for the Intensive Mode, 


rJ(0) = K2(2ir0/)Pp|(T+f-t^) 


2 2 
0 = 80/ 
1 


The impulse response of the single pole RC video filter is 

xn 


( 46 ) 



( 47 ) 


h (t) = 

V 




The double convolution of (47) with (44) may .be acconq>lished in closed 
form with the foiloxring result; 


R^(T)eh^(T)fih^(T) 


R^(0) 


~1 


1 

2 



2 \ 


2 2 


exp^a2T+a2''cr 





(48) 


Equation (48) may be rexnritten in a more simple form; 

rs-r * 

2 

where for Global Mode Cy=0.842 and for Intensive Mode ?^-0. 783 and P^(x) 
is normalized to a maximum value of one at T = 0. 


Rexjriting the noise and signal autocorrelation functions in (39) in 
a form similar to (49) , the variance of the Integrating gate output is 
given by 




[‘ 



dT 


P^(T)P^(T) 


(50) 
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where, from equations (41) and (43) , = 0* ^18 for the Global Mode and 

r =0. 964 for the Intensive Mode. The noise autocorrelation coefficients 
n 

in (50) are equal to the curly. bracketed factors in equations (40) through 
(43) . Figures 8 through 13 compare some of the autocorrelation coefficients 
in (50) for both before and after filtering. Figure 8 shows, for the Inten- 
sive Mode, that the video filter has little effect upon the noise alone and 
since the signal decorrelates more slowly than the noise, the assmaption in 
(37) is justified. Figure 11 illustrates that for the Global Mode there is 
a significant difference in the noise and signal decorrelation times. Fig- 
ure 8 also indicates the degree to which the radar in the Intensive Mode 
departs from an ideal matched filter, i.e. for matched filter conditions, 
p^(x) = particular note In Figures 8 through 13 is the fact that 

the video filter has a much more dramatic impact upon the noise than on the 
signal, especially, in Global Mode. 

2.3.5 Definition of Signal-To-Noise Ratio (SNR) 

2 2 

In Equation (50), it is tempting to call the ratio R^(0)/R.^(0) the 

square of the signal- to-noise ratio; however, this point deserves some 

2 

discussion. First, it must be remembered that R^ (0) is not the same for 

the Plateau and Attitude/ Specular gates; that is, the square of the aver- 

2 

age return power, i.e. R.^ (0) , evaluated at the center of the Plateau gate 
is certainly different from the average power at the midpoint of the Atti- 
tude/Specular gate. This fact is illustrated in Figure 14 which shows how 
the ratio varies as a function of pointing angle and mode. This, figure 
merely demonstrates the fact that as the pointing angle increases the peak 
in the average return power occurs later in time due to the increased delay 
1 -^Tnp along the antenna boresight axis to the surface. In fact, as the 
pointing angle approaches or exceeds one beamwidth, the time at which the 
maximum occurs is given by 


t 



This point raises the interesting question of how the altimeter signal-to- 
noise ratio should be defined! Usually, the term "signal" is taken to be 
the peak of the average return power; however, such a definition can be 
very misleading when considering large pointing angles. For example, when 
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CoBjpaxlson of the autocorrelation coefficients 
of the signal (no video filtering) and noise 
(before and after video filtering) for the 
Intensive Mode. 
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Figure 9. Comparison of the autocorrelation coefficients 

squared of noise before and after video filtering 
for the Intensive Mode. 
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NORMALIZED AMPLITUDE 



Figure 10. Comparison of the autocorrelation coefficients 
squared of the signal before and after video * 
filtering for the Intensive Mode. 



TIME T (IN ns.) 


Figure 11. Comparison of the autocorrelation coefficients 
of the signal (no video filtering) and noise 
(before and after video filtering) for the 
Global Mode. 
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Figure 12. Comparison of the autocorrelation coefficients 
squared of the noise before and after video 
filtering for the Global Mode. 
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Figure 13. Conq)arison of the autocorrelation coefficients 
squared of the signal before and after video 
filtering for the Global Mode. 
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figure 1‘4. The ratio of the Plateau gate power to the Attitude/Specular 


2.4 


the pointing angle is small, the Plateau gate will measure the peak of the 
average return power. As the pointing angle increases, the Attitude/Specu- 
lar gate will measure the peak and, finally, as the angle increases st ill 
further neither the Plateau gate nor Attitude/ Specular gate will encompass 
the peak of the average return. In the latter case, the peak return power 
is not meaningful since it does not tell what the power level is in the 
vicinity of the tracking gates. Since maintaining track-lock is the most 
important fiinction of the altimeter, it is the power in the Plateau (and 
Ramp) gate which is most relevent. For this reason, the term "signal" (as 
used in signal-to-noise ratio) will be defined as the average return power 
at the midpoint of the Plateau gate. 


2.3.6 The Error for a Ten Second Average 


With the above definition in hand, equation (50) should be rewritten 

as two equations; one for the Plateau gate and one for the Attitude/Specu- 

2 2 

lar gate. For the Plateau gate, define R^(0) = R^(0), and, similarly, 
for the Attitude/ Specular gate, define r2 (o) = R^(0), then with 

t 2 2 "1 1/2 ^ 

R „(0)/R (0) I , equation (50) can be rewritten for the Plateau 
vp n j 

and Attitude/Specular gates as follows; 
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For high SNR ( ~ 20 dB) , the different power levels in the two gates has no 
appreciable impact on the variances of the gate outputs. However, for low 
SNR and large pointing angles the difference in gate power levels tends to 
suppress the variance of the Attitude/ Specular gate output as well it should 
since the Attitude/Specular gate "sees" the higher power level (see Figure 
14). 

Equations (51) and (52) determine the variance of the output of the 
integrating gates. The variance of the output of the one second averager 
along with the additional ground based averaging (to form a ten second aver- 
age) is given by 


Var(e „ 
as 


) ~ 


10 M 
as 


(53) 


and 


Var(e ) z 

IT 


Var(e ) 
_E_ 

lOM 

P 


(54) 


Using the same technique as in Section 2.2 to determine the average estlma-' 
tion function, the variance of the estimation function is fotmd to be as 
follows; 


Var(A) - 
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(55) 


or, equivalently. 


Var(A) = 


e /G 
as as 


e /G 
P P 


2 r 


Var(e ) 
as 


Var(e ) 


lOM (e )^ ICM (e )^ ^ 
as as p p 


(56) 


The average integrating gate output is given by equation (8), however to the 
same order of approximation used to find the variance of e, (8) can be ap- 
proximated by the following; 
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( 57 ) 


e » GR (0) 

V 

Substituting (51), (52) and (57) In (56) yields a result for the variance 
of the estimation fxmction which involves three easily performed numerical 
integrations for each mode. 

Using the curves in Figure 6, the variance estimates for A can be 
translated into variance estimates for the pointing angle ^ for a ten sec- 
ond average. Figure 15 shows the error as a fxmction of pointing angle and 
SNR while Figure 16 shows similar results for the Intensive Mode. In com- 
paring these two results, it is very interesting to note that although the 
Global Mode Plateau gate averages sixteen times as many returns per xmit 
time as does the Intensive Mode, the resultant error for the two modes and 
high SNR are comparable. This is due to the fact that the estimation func- 
tion curve for the Intensive Mode has a much greater slope than the Global 
Mode (see Figure 8). These curves clearly show that this technique for 
estimating the radar altimeter pointing angle is a most powerful approach. 

2.4 Other Error Soixrces 

The previous section has been concerned with the errors which are due 
to the random nature of the backscattered return and the Influence of noise. 
In addition, there are two other sources which can significantly contribute 
to an erroneous estimate of the pointing angle. 

2.4.1 a® Variation With Angle of Incidence 

The first of these other error sources involves the angular behavior 
of CJ® for very near normal incidence. Throughout the previous analysis it 
has been implicitly assumed that O® is constant over the range from zero 
degrees to the equivalent angle of the Plateau and Attitude/ Specxilar gates. 
That is, in the equation for the flat surface impulse response, i.e. equa- 
tion (10), Cf®(x|»^) has been assumed to be constant. Table III shows the rel- 
ative time of occtxrrence of the midpoints of the gates and the corresponding 
angles of incidence on the mean flat surface. Since the greatest angle is 
about one degree, the above assumption on O® is certainly not unreasonable. 
Previous radar altimeter measurements of a a® [7, Chapter 8] indicate that 
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Figure 15. The precision of the pointing angle estimation 
technique for Global Mode. 
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Figure 16. The precision of the pointing angle estimation 
technique for Intensive Mode. 


TABLE III 


The Angle of Incidence Correspondiag to the Midpoints 
of the Plateau and Attitude/ Specular Gates 



INTENSIVE MODE 

GLOBAL MODE 

Plateau 

Gate 

Attitude/ 
Specular Gate 

Plateau 

Gate 

Attitude/ 
Specular Gate 

Mid Point 
Relative Time 
(ns) 

68.75 

800 

400 

800 

Corresponding Angle 
of Incidence (deg) 

0.28 

0.97 

0.68 

0.97 


for the most ocean surface conditions, the assunq?tion of a constant <J“., 
out to about one degree, is valid. Scattering data obtained by GEOS-3 
have indicated that there are, however, cases where this assui^tion is 
violated. Such situations . have been shown to occur when the surface has 
a very gmall mean square slope [7, Chapter 9]- Since a small mean square 
slope generally implies a lack of high frequency components in the surface 
height spectrum, ocean surfaces giving rise to this type scattering are 
possibly either void of any significant wind field or are covered by a 
more viscous material such as oil. 

If a“ decreases over the equivalent angle of incidence range represent- 
ed by the integrating gates, this can lead to an underestimate of the point- 
ing angle. That is, if a° decreases with angle of incidence, the trail- 
ing e4ge of the average return will decay faster than if 0° were constant. 
This increased decay, however, could very well he interpreted as a smaller 
pointing angle error. Of course, when the roll-off of O® with angle of 
incidence is sufficiently rapid as to produce a return whose trailing edge 
is below the nadir return (for 0® assumed constant), the computed estima- 
tion ftinction will exceed the maximum possible value given in Figure 6. 

This situation has been observed in the GEOS-3 data and is easily spotted 















because the received power level Increases significantly and the altitude 
data noise level decreases. The more troublesome case where roll-off 
can result in an erroneous but nonzero estimate of the pointing angle ran 
be spotted by observing the change in the trend of the pointing angle data. 
That is, for such cases, there will be a distinct change in the behavior of 
the pointing angle estimates which is not consistent with spacecraft dynam- 
ics. It should be noted that the estimation technique is generally not ap- 
plicable over terrain due to the Inhomogeneous nature of the scattering 
surface. In conclusion, a rapid roll-off of a° with angle of incidence 
does not invalidate the technique but it does result in the possibility of 
erroneous estimates; however, suspicious data can be checked by methods 
discussed above. 

2.4.2 Altimeter Calibration 

All portions of the radar receiver modeled in this analysis are analog, 
see Figure 3, and are therefore subject to thermal and ageing effects. 

Since spacecraft power was not available to maintain a constant thermal en- 
vironment, extensive thermal/vacuum testing and calibration of the altim- 
eter were accomplished [6] . Because of latmch constraints, this calibration 
was not quite as extensive as would have been desired under a less severe 
schedule. In regard to calibration of the Plateau and Attitude/Specular 
gate outputs, two problems have come to light since launch. 

The first of these' two problems comprised inadequate thermal/vacuum 
calibration of the Global Hode Plateau Gate. For this gate, input /output 
calibrations were accomplished only at room ambient temperature and pres- 
sure. As a consequence, the data correction processing performed by the 
Wallops Flight Center essentially accomplishes a single point thermal cor-” 
rection for the output of the Global Mode Plateau gate. A cursory examina- 
tion of the temperature dependent calibration data in [6] for the Attitude/- 
SpectLLar and Noise gates* clearly indicates that thermal effects are cer- 
tainly not negligible. Unfortunately, there is no obvious means for over- 
coming this problem short of conducting special inflight tests in which the 


*The 200 ns Attitude/ Specular and Noise gates are common to both the Inten- 
sive and Global Modes and both were calibrated more extensively during. In- 
tensive Mode testing. 
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altimeter is switched back and forth between the Inteiisive and Global Modes. 
Assuming that the intensive Mode 'is operating within its thermal calibra- 
tion range, it should provide temperature independent estimates of the point- 
ing angle. By comparing pointing angle estimates obtained from both modes, 
it may be possible to generate a pointing angle correction curve for the 
Global Mode. Such testing is recommended when operational considerations 
permit such a configuration of the altimeter. 

As a sample of what effect inadequate calibration of the Global Mode 
Plateau gate can have, let the thermal environment be such that the effec- 
tive gain (from detector output to one second averager output) changes by 
5% from its preflight value. The average estimation ftanction is given by 
(21), i.e. 


A = 1 


e /g 
as as 


yx 

e /G 
P P 


(55) 


^ /s 

where G^^ and are the gains obtained from preflight testing. However, 
according to (8) the average gate outputs may be written as 


and 


e — G f 
as as as 


e 

P 


G f 
P P 


\ 


where G and G are the actual gains at the operating temperature of the 
as p ^ 

tracking gates. Taking the differential of (55) with G - G , i.e. no 

as as 

error for the Attitude/Specular gate, the change in the average estimation 
fxinction corresponding to a change in is given by 


dA = (1-A) 



(56) 


where if G -G the dG =0. It should be noted from (56) that if G is 
P ^ P P‘ _ P 

larger than then the change in A is positive and the pointing angle es- 
timate will be low, i.e. too small. Conversely, for a decrease in gain, 
the pointing angle estimate will be too large. Figure 17 illustrates the 
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POINTING ANGLE ERROR (DEGREES) 



Figure 17. .Pointing angle bias error for a 5% change in 
Plateau Gate gain for Global Mode. 


■pointing error resulting from a constant 5% change in the gain of the Global 
Mode Plateau gate.* If the gain increases by 5%, the computed estimation 
function will be greater than its maximum possible value of 0.306 for 
5<0.4“, and no estimate of the pointing error is possible. For either an 
increase or a decrease in the gain, the error in the computed pointing 
angle is always greatest near ^ = 0“ . Figure 18 is a plot .of estimated ver- 
sus true pointing angle for the 5% gain change. 

A preliminary comparison of Global Mode pointing angle estimates with 
spacecraft derived attitude data from [2] indicates that thermal calibra- 
tion may be a problem. That is, a plot of the radar estimated pointing 
angle versus the spacecraft sensor estimated angle is very similar in form 
to the upper curve in Figure 18. The mean error in the vicinity of zero 
pointing angle is higher, i.e. on the order of 0.7°. This would imply a 
gain increase of about 15%; however, the problem requires much further 
study. In particxilar, the data will be grouped by temperature to determine 
if, in fact, the temperature is a meaningful parameter in terms of data re- 
peatability. In any case, there is a strong indication in the Global Mode 
pointing angle data reduced so far that inadequate gate calibration may be 
a major source of error. 

The second calibration problem is similar to . the first but it involves 
the Intensive Mode and is not as severe as the first in regard to its im- 
pact. During the very early operation of the Intensive Mode, the altimeter 
tended to run somewhat colder than was expected. Consequently,, data were 
obtained at temperatures which were below the preflight calibration range 
and, thus, were only approximately corrected. Fortunately this problem will 
probably only result in the loss of some very early data since the operat- 
ing temperature has increased back into the calibration range. Furthermore, 
since it is known when the data are outside of calibration range, i.e. by 
comparing the operating temperature with temperatures during calibration, 
one knows exactly when to suspect the data. Preliminary comparisons of 
Intensive Mode pointing angle estimates with spacecraft sensor derived at- 
titude data indicates a much closer agreement than in the case of the Glob- 
al Mode. 

*The pointing angle e^ror is defined- as the difference between ^he ttue 
angle (computed from A) and the incorrect an g le (computed from A + dA). 
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ESTIMATED POINTING ANGLE (DEGREES) 



TRUE POINTING ANGLE (DEGREES) 


Figure 18. Estimated versus true pointing angle for a 5% 
change in Plateau gate gain in Global Mode. 
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3.0 the impact of pointing error 

* The previous sections of this report have presented a technique for 
computing the radar antenna pointing angle and have estimated the accuracy 
and precision of the method. Given the fact that problems associated with 
inadequate calibration can be minimized, as per the discussion in the pre- 
vious section, the method has the potential of providing very accurate 
pointing data. The need for such pointing estimates, for this mission, 
rests solely upon its potential to increase the accuracy of the primary 
altimeter data such as altitude, surface scattering cross section and wave- 
height estimates. Since the GEOS-3 altimeter antenna beamwidth is relative- 
ly large, the need for pointing angle correction is not as great as in the 
case of, say, Skylab [7]. However, this does not imply that it is a negli- 
gible effect; the corrections for pointing angle effects depend entirely 
upon the attitude excursions encountered during inflight operation of the 
altimeter. This section will consider the effects of pointing errors upon 
the following quantities; the normalized surface scattering cross section 
or a®, the altitude data, and estimates of the ocean surface waveheight. 

3.1 cr® Estimation 

For the Global Mode, the average return power at the output of the ^ 
radar antenna is given by 

CO 

_ . I 

P (t) = T exp (- ^ sin I P (t-t’)a®(i|i’) 

4(4ir)"^L^Lph^ ' Y ^ 

0 

• [• ^ (y /? 

where P (*) is the transmitted waveform and' 

T 

tani|; = j^ct’/hj (58) 

Since P (t) is a function of relative delay time, the question arises as 
r 

to exactly where (in t) the measurement is made* Since the AGC control volt- 
age, during tracking, is a filtered version of the Plateau gate output [6], 
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the received power measurement < will be given by the following; 

500 




y'p (t) exp[(500-t)/800] dt 
300 ^ 


( 59 ) 


500 


f exp[(500-t)/800] dt 


300 


where time has the units of nanoseconds (see Figure 2). 

That is, the actual power measurement is taken to be an exponentially 
weighted average of the received power over the Plateau gate. Assuming 
that a® ) is reasonably constant over the equivalent angular width of 
P^(t), it may be removed from inside the integral in (57), i.e. 


P^(t) » ; §- 

4(4ir) L L h 
r p 


uu 

4 


exp ^sin^5j I P^(t-t')exp - — t’cos2sj 


(f sln25)df 


1 1/2 


where tan ip = |^ct/h^ . Rewriting (60) in the following form; 


(60) 


P (t) = P„ 
r T 


J 


Fg(t,5) 


(61) 


the average received power is given by 

500. 


<P^(t)> « P^ 


L "p . 


f Pg(t,0 exp [(500-t)/80o] dt 


300 


500 


f ^ [(500-t)/8O0] dt 


300 


(62) 


where corresponds approximately to the midpoint of the Plateau gate, i.e. 
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z 0.68®. Since <^P^(t))> and are measured and ,a®(ip^) Is to be com- 
puted, a plot of the ratio 


<Pr(t)> 



is convenient because it presents a function of received 

power, transmitted power and pointing angle. Figure 19 illustrates how 
the ratio depends upon pointing angle and may be used in conjunction with 
measurements of P (t) > and P„ to determine fcf®(^ )/L J . The curves in 
Figure 19 were generated using the following parameters obtained from [6]; 

G * 36.5 dB 
o 

X *= 0.02158 m 
c = 0,3 m/ns 

L = 0.8 dB 
r 

h « 843 km 

-3 

Y = 1.49x10 (2.4® antenna beamwidth) 

It should be noted that the results for the Global Mode in Figure 19 should 
not be used when there is an indication that 0® (up) may be changing rapidly 
with incidence angle. 

For the Intensive Mode, the development follows that given for the 
Global Mode with exception that the waveform, P^(*), requires more careful 
definition. The actual transmitted waveform comprises a 1.2 ps pulse 
whose carrier frequency is linearly swept at a fixed rate. The post-IF 
portion of the radar receiver operates on a compressed or deramped version 
of the expanded average return. Thus at first glance and as far as the 
measurement of 0® is concerned, it would appear that F^(*) in (57) could 
be replaced by the ideal post compression waveform. With fhe exception of 
time sidelobes, this statement is true! That is, after compression in the 
receiver, the actual predetection point target response has, essentially, 
no time sidelobes, see Figure 4. Thus, for an extended target, the normal 
build-up of- power due to integrated time sidelobes does not occur and the 
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Plateau gate responds only to the Integrated mainlobe of the point target 

response. Since the Plateau gate provides the AGC control voltage, from 

which the received power is inferred, the time sidelobes should be ignored. 

Therefore, for the Intensive Mode, P^(t-t’) in (57) can, equivalently, be 

replaced by 100 •Pjj.(t-t’) where 100 is the ideal compression gain* and 

p (t-t’) is the IF point target response. This approach is different from 
IF 

the scheme used to calibrate the AGC but is entirely equivalent** to it, 
see [6]. 


With 


gVc(IOO) 

4(4ir)^L^h^ 




(t-t') exp 


0 


[-||tVos25] 




(63) 


the mean return power, time-averaged over the Plateau gate, is given by the 
following ; 


<P^(t)> 



75 

/ F^(t,5) exp [(75-t)/50] dt 
62.5 
75 

f e^[( 75-t)/5o] dt 
62.5 


(64) 


where is approximately equal to 0.28®* Figure 19 shows how the ratio 


<p (t)> 



varies with pointing angle for the Intensive Mode; 


*The AGC calibration accounts for any differences between ideal and actual • 
compression gain. 

**The difference occurs because ideal pulse compression devices do not 
physically exist. 
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■The curves in Figure 19 essentially illustrate the degree of correc- 
tion that is necessary to account for the effect of pointing error in com- 
puting a° from the received power data. If the pointing error is less than 
0.8”, the maximum correction for the Global Mode is 1.75 dB while for the 
Intensive Mode it is 2.14 dB. 

3. 2 Global Mode Altitude Bias 

The GEOS-3 altimeter employs a so-called split gate tracker to locate 
and follow the backscattered return. For such a tracker, the time of oc- 
currence of the Ramp and Plateau integrating gates is constantly, adjusted 
until, on the average, the following tracking law is satisfied; 

2e - e =0 (65) 

As shown in Figure 2, for an idealized return with a linear rise equal to 
6 and the gate widths also equal to 6 and no pointing error, the leading 
edge of the Ramp gate will occur 2h/c seconds after detection of the tr^s- 
mitted pulse. As the pointing angle increases, the trailing edge of the 
return departs from its ideal behavior which causes e^ to change and- the 
tracking point must also change. Thus, the leading edge of the ramp gate 
no longer occurs at 2h/c, and the altitude estimate provided by the altim- 
eter is biased. A bias will also result if the Ramp gate is not matched to 
the rise time portion of the average return waveform- 

Figure 20 shows that the altitude bias for Intensive Mode is not a 
very sensitive function of pointing angle. This results from the fact 
that the Ramp gate width (12.5 ns) is very well matched to the average 
return rise time and the Plateau gate is reasonably insensitive to point- 
ing errors. For the Global Mode, the situation is considerably different. 
Figure 21 illustrates the normalized video return presented to the track- 
ing gates during the Global Mode for nadir pointing. Of particular note 

t 

is the rather long rise time of the leading edge which is due, in great 
part, to the long trailing edge of the video point target response (see 
Figure 5) . The nominal midpoint of the Plateau gate for an idealized re- 
turn ( - 400 ns) occurs when there is a rather significant droop in the 
trailing edge of the return. Therefore, it is obvious that the Global Mode 
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Figure 20. Intensive Mode waveform induced bias as a 

function of pointing angle (for H, ... < 0.5 m) 



bias will be significantly different from the Intensive Mode. 

Figure 22 shows the altitude bias as a function of pointing angle for 
the Global Mode. According to the definition of bias as given in the fig- 
ure, the altimeter will always produce an altitude measurement which is 
too large, i.e. relative to 2h/c. Although not noted explicitly in the- 
derivation of the estimation function and its variance, this shift in the 
-location of the gates was properly accounted for in the actual computa- 
tions. The results in Figure 22 clearly show that bias for the Global 
Mode is not negligible. It should be noted that the bias results present- 
ed in Figure 22 are only for the waveform and pointing induced effects; 
Reference [6, page 203] presents other bias factors which result from the 
GEOS-3 design. Reference 6 has also computed the waveform induced bias at 
nadir and obtains a value of 4.95 m as opposed to the 6.15 m shown in Fig- 
ure 22. The difference is attributed to a slighly different result for the 
video average return waveform, i.e. compare Figure 20 with Figure 85 (page 
205) of Reference 6. 

3.3 Waveheight Estimation 

According to the theory of linear random surface scattering, surface 
roughness tends to increase the rise time of the average return because 
of the range distributed nature of the specular points on the surface. To 
avoid the mathematical horrors associated with numerical deconvolution, one 
is usually content with measuring the increased rise time of the return 
and then translating this into- an equivalent rms surface height via some 
straightforward model. The model normally requires an accurate computation 
of the so-called flat surface response, i.e the rms surface height is very 
small relative to the range extent of the video point target response. For- 
tunately, in the case of GEOS-3, the effect of pointing error upon the rise 
time extent of the flat surface response is negligible until the pointing 
angle approaches 1.2® (one half a bearawidth) . For this reason, pointing 
errors can usually be ignored in waveheight estimation. However, if the 
pointing error does approach 1.2®, it should be properly accounted for or 
it will give rise to an overestimate of the surface waveheight. That is, 
pointing error tends to increase the rise time of flat surface return. 
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POINTING ANGLE (DEGREES) 


Figure 22. Global Mode waveform induced altitude bias 
as a function of pointing angle. 
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4.0 AGC CALIBRATION FOR FLUCTUATING TARGETS 

As noted in section 3.1, the determination of a® assumes an accurate 
measure of the average return power at the Plateau gate for both the Inten- 
sive and Global Modes. For the GEOS-3 altimeter, this measurement is ac- 
complished by sampling a filtered replica of the AGC control voltage and, 
using preflight calibration curves. Inferring the return power. The ac- 
curacy of this approach is, therefore, highly dependent upon the accuracy 
of the preflight calibration curves. 

4.1 The Clean Versus Clutter Problem 

Preflight calibration curves of AGC voltage versus average received 
power were obtained over a wide range of receiver temperatures and received 
power levels, in all cases, the average modulation waveform was tailored 
to simulate the expected inflight average return waveform. In addition, 
data were also obtained for the case where the basic waveform was, essen- 
tially, ftirther modulated by a random noise source so as to simulate the 
fading and fluctuating nature of the backscattered return. When the simu- 
lated average return was not noise modulated, the data were referred to as 
the Clean AGC calibration curves; when the simulated retturn was noise modu- 
lated, the data were called Clutter AGC calibration curves. Although not 
representative of actual backscattered retijms, the non-noise modulated or 
Clean AGC calibration data were obtained to cross check measurement techni- 
ques and basic AGC stability. 

Since the AGC control voltage was determined by the power in the Pla- 
teau gate, it seemed reasonable to compare the Clean and Clutter AGC data. 
Because the Clutter data were obtained by averaging a large number of re- 
turns, the statistical error in the Clutter data was expected to be small; 
thus, the Clean and Clutter data were expected to be nearly identical. 
Figure 23 shows a comparison between Global Mode Clean and Clutter AGC data 
for nearly identical receiver temperatures. The agreement of the data is 
very good, the average difference is about 0.4 dB and this is certainly 
within measurement tolerances.* However, the comparison for the Intensive 


*The majority of this 0.4 dB difference is due to the 0.6®C temperature 
difference between the two sets of data. 
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(volts) 



Figure 23. Global Mode Clean and Clutter AGC 
calibration data, from [6]. 
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Mode was quite poor as shown in Figure 24. For these data, the differ- 
ences ranged between 3.5 and 4 dB with the Clean data always indicating the 
higher voltage for equal Glean and Clutter power. 

Extensive verification of the test procedures and techniques for the 
Intensive Mode Clutter data was accomplished by E. L. Hofmeister and B. N. 
Keeney of GE, but no obvious problems were found. The disparity was par- 
ticularly troublesome because a decision had to be made as to which set of 
data for the Intensive Mode should be used to convert inflight AGC voltage 
measurements into received power. The Clean AGC curves were chosen to re- 
duce the inflight date since the basic calibration process was considered 
to be more accurate. However, a preliminary analysis of cr“ data resulting 
from both the Intensive and Global Modes showed that the Intensive Mode, 
using the Clean AGC curves, was yielding estimates of 0® which were con- 
sistently 3 to 4 dB below the estimates obtained from Global Mode. If, 
however, the Intensive Mode AGC voltage data were converted to received 
power using the Clutter curves, the 3 to 4 dB difference could be reduced 
to less than 1 dB. Thus, inflight data seemed to indicate that the Clutter 
calibration data should have been used to convert AGC voltage to received 
power . 

Although the comparison of inflight a® values resulting from both 
altimeter modes is a reasonable means for resolving the Clean vs. Clutter 
question, it still left much to be desired. In the first place, the error 
bounds in such a comparison mught be comparable to the 3.5 to 4 dB differ- 
ence which was to be resolved. That is, because of varying surface condi- 
tions, spacecraft attitude and other factors, comparisons of this type are 
always subject to uncertainties. In addition, even if 0® comparisons did 
show which set of calibration data should be used, it would not answer the 
basic question of why the difference in the first place. -Because future 
altimeters will also probably employ AGC systems of the GEOS type, it seem- 
ed imperative to find the source of the discrepancy. Since the calibration 
process was thoroughly checked and verified, it appeared that the discrep- 
ancy might be due to some fundamental difference between how the altimeter 
responded to a deterministic and a random signal. For this reason, the 
basic operation of the AGC loop was more thoroughly investigated. 
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Figure 24. Intensive Mode Clean and Clutter AGC 
calibration data, from [6]. 
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4.2 A Self-Consistent Solution To The Clean vs. Clutter Problem 

A solution to the Clean vs. Clutter problem must satisfy the follow- 
ing conditions. First, the solution must comprise no unjustified assump- 
tions in regard to the physics of the process. Second, as the experimen- 
tal results in Figure 23. indicate, the Clean and Clutter AGC curves for 
the Global Mode must nearly overlap. Third, the solution must very nearly 
account for the discrepancy between the Clean and Clutter data for the In- 
tensive Mode. It should be noted however that, in the final analysis, the 
proposed solution must be verified by tests and measurements on the actual 
hardware ! 

The results in Figures 23 and 24 indicate that for a non— fluctuating 
received power (Clean case), the AGC control voltage V may be expressed as 
a logarithmic function of the received power,* P^, i.e. 

V 2 a + b [lOlogCPp] (66) 


Equation (66) is strictly true only over a limited range of P^ values; how- 
ever, for the present, it serves to illustrate some very important points. 

It is necessary, next, to know how V responds to a fluctuating received 
power on a sample-by-sample basis rather than on the average as presented 
in Figures 23 and 24. To a first order, at least, (66) should nearly rep- 
resent the sample -by-sample case but with V dependent upon some finite pulse 
average of the received power (due to filtering in the AGC loop) . Although 
it is more reasonable to represent this average as a weighted average of 
the pulse-by-pulse power, the problem will be simplified to an equivalent 
uniform average of discrete values. This assumption will not significantly 
alter the results to be presented. The next question is as follows; should 
(66) be replaced by a "finite pulse average of logs", i.e. 


= a + b 


m .. 


or a "log of the finite pulse average", i.e. 


( 67 ) 


*The symbol log will be used for log- „ while £n will be used for log . 

XU G 
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( 68 ) 



Using the Clean AGC data there is no way to resolve this question and the 
AGC details in [6] are also insufficient in this regard. Both of the above 
equations merely state that any observation of the AGC control voltage is 
based upon < a finite pulse average of received power. 

Consider for the moment that either equation (or both) may be true. 

First of all, it must be realized that- although was explicitly written 

in both equations, this is incorrect for a sample-by-sample description of 

the control voltage. That is, the AGC loop is not excited by the received 

power, per se, but rather by a voltage which is proportional to the output 

of the Plateau Integrating gate, see Figure 25. Since the Plateau gate is 

nearly "matched" to the video point target response, the statistics of the 

gate output are virtually identical to the statistics of input voltage.* ** 

Thus, the output of the Plateau gate is nearly equal to (GP /K) where K is 

r 

the AGC attenuation, G is the system gain from antenna output to Plateau 
gate output (excluding the AGC attenuator), the probability density of P 

r 

is exponential, i. e. 


f(P ) =-i-exp(--^)H(P ) , (69) 

\ KP^/ 

and is the average received power at roughly the mid point of the Plateau 

gate. For the Intensive Mode, GP^/K is fed directly to the AGC loop for 

each ptilse. For the Global Mode, sixteen returns resulting from a single 

pulse burst are first averaged and then fed to the AGC loop. Thus, the 

AGC block diagram in Figure 25 may be replaced by its equivalent in- Figure 

26. Figure 26 clearly shows that the effective input to the AGC loop is 

GP where for the Intensive Mode P » P , and for the Global Mode 
r r r 


*This statement is theoretically demonstrated, for the first two moments, 
in Section 2 for both altimeter modes, and experimental verification for 
the Intensive Mode is given in [6, page 145]. 

**This fact is ignored in Section 5.2 (AGC Analysis) of [6], however it is 
clearly stated on page 90 of [6] in text and is shown in the second sheet 
of Figure 7 of [6] . 


63 




TO TRACKING 
GATES 


Figure 25. Block diagram of the AGC loop. 





Also, ia equations (67) and (68), P^(i) should now be replaced by P^(i-), 

The reason for replacing the true AGC system by the equivalent system 
in Figure 26 is to emphasize the fact that in the case of the Global Mode 
the AGC input is actually a sixteen pulse average while the 3iitensive Mode 
input is not averaged. Thus, the probability density ftmction for is 
exponential for the Intensive Mode and is given by the gaTnrua density [14] 
for the Global Mode,-i.e. 


f(P ) = . 
r 


^ 

y roj) 


(Pr)’^^e::p(-HP^/P^)D(p^) 


(71) 


where N = 16. This sixteen pulse averaging is the reason why the statisti- 
cal fluctuations in V are much lower for Global Mode than for the Inten- 
sive Mode. 


Using the gamma density for the Global Mode the average of the AGC con- 
trol voltage can be computed for both possible forms of V, i.e. equations 
(67) and (68). Since the average of the control voltage for the Clutter 
case must equal the Clean results, both (67) and (68) must satisfy the fol- 
lowing identity for the Global Mode; 


V = a 



(72) 


For (67), the average control voltage is given by 



Sxnce the P^(x), x 1,2, • * »,m, are independent 


*The double indexing of for the Global Mode is necessary because of the 
pulse burst operation, in (70) , i denotes the burst number while j denotes 
the pulse number within a burst- 
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where k = 4.3429. The probability density ftmction o£ y = JJn(P^) is easily 
shown to be 


a + 10b E 
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Thus, the mean value of V is given by 


V^ = a + 


_^kbH^ f 

(p.)”r(N) J 


y ejcp 


N 

Ny - — exp(y) 


dy 


(74) 


Integrals of this form are evaluated in the ^pendix with the following 
result; 


V = a i — 

1 ^ ^ r(N) 


[r(N)£n(P^) - r(N)iln(H) + T’ (N)] 


(75) 


or 


= a + kb[£n(P^) - fei(N) + 

Equation (76) may be simplified to the following; 

= a + b 1 10 log (Fp - k£n(N) + kifiCN) 


(76) 


(77) 


where if((N) is the digamma function which is defined for integer values of 
K as follows; 


ip(N) = 


- y 

N-l . 
k=l 


N = 1 
N > 2 


and Y is Euler’s constant, i.e. Y = 0.5772. With P^iag(N) = ~ k«.n(N) + ki{i(N), 
(77) laay be written as 
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V^ = a+b{l01og(pp 


(78) 


Equation (78) Indicates that the average AGC control loop voltage dif- 
fers from the deterministic or Clean case by the term bP^^ .^^ (N) . Stated 
another way, (78) implies that if the Clean curves were used to estimate 
the average fluctuating power, one would obtain a biased est ima te of the 
average received power. This statement merely reflects the well known fact 
that a logarithmic system (such as the AGC loop) produces a biased esti- 
mate of the average fluctuating power [15]. The fact that the bias is a 
function of the filtering prior to the AGC loop is shown in Table IV . For 
N = 16, the bias is -0.14 dB; this means that if eqviation (67) is the proper 
description of the AGC loop control voltage, then the Clean and Clutter data 
in Figure 23 shotild differ by no more than 0.14 dB. The data in Figure 23 
clearly show this degree of agreement. Thus, equation (67) produces a re- 
sult which agrees with both the Clean data and the average Clutter data 
for the Global Mode. 


TABLE IV 

AGC Loop Bias As A Function Of The Number Of Independent 
Returns Averaged Prior To The AGC Loop 


N 


1 

-2.51 

2 

-1.18 

3 

-0.76 

4 

-0.57 

8 

-0.28 

12 

-0.18 

16 

-0.14 

20 

-0.11 
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If the control voltage is as described in (68) , the average is given 


by 


or 



?,(«]} 




P 

r 



Since the P^(j»i) are independent for each i and j, the probability densi- 
ty of X, where 



j=i 


is given by 


f(x) = 


(Nm) 


Nm 


(P^)”“r(nM) 


(P^)^”l exp [- (HinP^)/?^] U(P^) 


(79) 


Thus, the average control loop voltage is given by 

Vj - a + b {lO log^oCP^) + 

where the analysis for has been used to deduce this result. Since 
ra > 10 for the Global Mode [6], the product Nm > 160 and neg- 

ligible. It is interesting to note that whereas the bias in V- depended 
only upon pre-loop filtering, i.e. N, the bias in V^- depends on both pre- 
loop and loop filtering, i.e. Nm. 

The results for and V 2 both agree, to within experimental error, 

with the data in Figure 23 for the Global Mode; that is, the Clean and Clut- 
ter curves are nearly identical. Thus, the problem of choosing which func- 
tion is a proper representation for the AGC loop control voltage still per- 
sists. Referring to Figure 25, the AGC control voltage V is a scaled and 
offset replica of the output from the loop filter. Since the filter is 
linear, the output represents the average of a number of inputs. If the 
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inputs are 10 log[P^(i)], i=l,2,...,m, then equation (67) clearly repre- 
sents a linear process which is characteristic of an averaging filter. Con- 
versely, there does not appear to be any obvious way for (68) to be the re- 
sult of a linear averaging process. That is, let q(P^(i)) be the ijA input 
to the loop filter in Figure 25. The output of the loop is then equal to 


hence, for (68) to be valid, the following equation must be satisfied; 

m 

= (81) 
1=1 

This equation has no solution for q^P^(i^ other than the trivial one which 
Implies no averaging by the loop filter. Thus,, based on the logarithmic 
behavior of the control voltage and the linear nature of the loop filter, 
it is concluded that is a realizable representation for the control volt- 
age and is not! For this reason, will be dropped from further con- 
sideration. 

For the Intensive Mode, the average AGC control voltage is given by 
the Global Mode result* with N = 1; that is, for the Intensive Mode, the. 
Plateau gate output is fed directly to the AGC loop on a pulse-by-pulse 
basis. • Thus, 

“ a + b|l01og(pp - 2.5l| (82) 

since the Plateau gate output is not filtered or averaged prior to its in- 
put to the AGC loop. Since the slope constant b is positive, (82) implies 
that the Clutter curves will be below the Clean curves as in Figure 24. 

The results in Figure 24 indicate that the difference between the Clean and 
Clutter curves varies from about 3.5 to 4 dB while (82) shows that the theo- 
retical difference should be no more than 2.51 dB. Thus, there still is 


10 log 


IE p^(i)l 


*The constants a and b are the same as may be verified by comparing IM and 
GM Clean curves. 
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some discrepancy which cannot be accounted for. With the thought that the 
residual difference may be due to the linear approximation in (66) , the 
control voltage was assximed to be more accurately represented by a quadrat- 
ic expression of the form 

+ b^l01og(Pp + Cj^[l01og(P^)J^ (83) 

for the Clean or deterministic data and 

= a 2 + b^ 10 log (P^) + c^ [lO log(P^)] ^ (84) 

for the Clutter or random data. The coefficients (a^,b^,c^) and (a^jb^jC^) 
were determined by a least square fit of the data in Figure 24 to the ex- 
pressions in (83) and (84) . Substituting 

igl01og[p^(i>] 

for 10 log (P^) in (84), the average was analytically determined. Compar- 
ing, the difference between (83) and (84) with the difference between (83) 
and V^, resulted in a voltage difference which translated into the same 1 
to 1.5 dB discrepancy as before. This fact merely confirmed the suspicion 
that the probability density of 101og^P^(i)J was so concentrated about 
10 log (P^) that a linear approximation was sufficient, at least on a piece- 
wise basis. 

In conclusion, 2.5 dB of the 3.5 to 4 dB difference between the Inten- 
sive Mode Clean and Clutter AGC calibration curves has been accounted for. 
This 2.5 dB difference results from the logarithmic relation between re- 
ceived power and the AGC control voltage and the lack of any filtering or 
averaging between the Plateau gate output and the AGC loop input. The 
difference does not exist in Global Mode because of the sixteen pulse aver- 
aging between the Plateau gate output and the AGC loop input. The analysis 
presented here clearly shows that the Clutter data should be used to con- 
vert inflight AGC voltages into received power, except for BIT/CAL where 
the Clean curves should be used. The residtial discrepancy of 1 to 1.5 dB 
cannot be accounted for at this time. If it could be attributed to 
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measurement error, the Clutter curves could be moved closer to the Clean 
curves by 1 to 1.5 dB to account for the error. However, given the com- 
plexity of simulating a noise-like coded rettim such as in the case of the 
Intensive Mode Clutter data, it is not completely clear what the source of 
the error is. In the author's opinion, achieving an accuracy of 1 to 1.5 
dB for such a system is a great credit to the engineers who designed, built 
and tested the GEOS-3 radar altimeter. It is anticipated that during the 
lifetime of the mission many opportunities will arise in which similar In- 
tensive and Global Mode data sets can be compared and analyzed with the 
purpose of resolving the final 1 to 1.5 dB discrepancy in the Intensive 
Mode Clutter calibration data. 

5.0 PROCESSING FOR AND INITIAL ESTIMATES OF POINTING ANGLE 

The altimeter antenna boresight angle with respect to nadir is a mat- 
ter of interest in the refinement of certain altimeter derived target fea- 
tures such as radar cross section and ocean waveheight as well as the alti- 
tude data. Though the degree of effect upon the accuracy of the various 
quantities of interest is variable, it is certainly not always negligible. 
The theoretical development of the estimation technique and its error bounds 
. is given in section 2 . This section will be concerned with the application 
of the formal restilts to actual GEOS-3 data. 

5.1 Altimeter Pointing Angle Estimation Algorithm 

Basically, there are only two measurements -which the GEOS-3 altimeter 
produces that are used in pointing angle computation; the one-second aver- 
aged Plateau integrating gate output and the Attitude/Specular integrating 
gate output. From these, the "average estimation function" is computed 
from the following; 


A = 1 -3^ (85) 

APG ■ 


where APG represents an average of the Plateau gate outputs (processed to 
remove the effect of the gate's transfer function gain) and ASG represents 
the same for the Attitude /Specular gate. For the GEOS-3 prf, ASG and APG 
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must be based upon at least a ten second average for (85) to correspond to 
the estimation function given by (19). Using the theoretical curves of A 
as a function of pointing angle shown in Pigure 6, (85) can be used to es- 
timate the altimeter pointing angle. 

The GEOS-3 telemetry system operates in two different data rate modes 
in which the Plateau and Attitude/Specular gate outputs are sampled at 
slightly different rates [17]. However, in both cases, the telemetry data 
rate permits an oversampling of the outputs since the basic gate outputs 
represent one-second averaged values and the telemetry sampling is at least 
two times per second. Processing of GEOS-3 data is largely organized on a 
per frame basis so it is convenient to define intervals in terras of frames. 
A low data rate major frame period is 2.048102 seconds, so five of these 
frames adequately cover a ten second interval. High data rate frames are 
3*276964 seconds long, so three of these comprise very nearly ten seconds 
of data* 


Nimerical approximation algorithms were developed to facilitate the 
transformation from A to ^ (in degrees). For Global Mode, the expression 
is as follows; 


5 = [s. 


14848 - 10.2796 Jin (A + 1.9033)] 


1/2 


( 86 ) 


where ? is restricted to be less than 2.4“ . For the Intensive Mode, the 
approximate expression is 


C = [ 5.0935 And. 9976 -A) “ 2.04346] 


1/2 


(87) 


and 5 should be limited to less than about 2° . Table V compares pointing 
angles computed from (86) and (87) with those obtained from Figure 6 as a 
function of the estimation function, A. 

GEOS-3 altimeter data presently exists in various forms, but the eas- 
iest to manage and transport is magnetic tape. This necessitates a com- 
puter program to access the data, find values of interest, aiui make the 
necessary computations. Such a program (in Fortran) has been developed and 
is currently running on the Wallops Flight Center ECLIPSE computer. The 
program utilizes a subroutine written by Dr. G. S. Hayne of Applied Science 
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TABLE V 


Comparison Of Pointing Angles As Obtained From Figure 6 
With The Numerical Approximations Of Eqtiations 86 and 87 


GLOBAL MODE 


A 

5 : Figure 6 

C : Eqn. -86 

0.306 

0“ 

0.01“ 

0.2725 

0.4“ 

0.396 

0.1736 

0.8“ 

0.798“ ' 

0.0148 

0 

• 

H 

1.205“ 

-0.1943 

1.6“ 

1.625“ 

-0.4362 

2.0“ 

2.05“ 

-0.6795 

2.4“ 

2.46“ 


INTENSIVE MODE 


A 

5 i Figure 6 

5 : Eqn* 87 

0.504 

0“ 

0. 0003“ 

0.458 

0.4“ 

0.393“ 

0.308 

0.8“ 

0.793“ 

0.020 

1.2“ 

1.196“ 

-0.472 

1.6“ 

1.6004“ 

-1.278 

2.0“ 

2.0“ 
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Associates, Inc . , which reads into core one entire frame of either high 
or low rate data from external data files on tape or disk. These external 
data dumps must be in the GAP or ARC form. 

This program is versatile in that either Global or Intensive mode data 
(at a high or low data rate) can be used. A test is made during program 
execution to branch to the appropriate program segments for each case. All 
the quantities of interest are tested for value bounds, and, if within 
bounds, are used in the computations. If not, a message designating an 
out-or-range condition is printed and the program skips to the next con- 
secutive ten second interval. At present, there is no overlapping as in a 
sliding-window average and each value reported by the program represents a 
completely unique, ten second Interval. The time shown in the printout is 
the start time for the middle frame in each ten second inteirval. Plateau, 
Attitude/ Specular, AGC, altimeter status, transmitted power, estimated 
pointing angle (from A), and cr® averages are also printed. Since each sta- 
tus word represents an average, anything other than a whole number will im- 
ply a status change during the respective ten second averaging interval. 

If this occurs during automatic tracking, a loss of track is thereby indi- 
cated. Sometimes, very strong returns may be encountered. In such cases, 
the resultant values of A will exceed an upper bound and no estimate of 
pointing angle can be made. Ifhen this occurs, an easily distinguishable 
value of -9999.99 is assigned to the pointing angle and a° estimates. In 
such cases, however, all the other reported averages are given, and often 
a very high AGC or a loss of track may be noted. Figure 27 presents an 
abbreviated listing of an output for Global Mode Pass 184 by this program. 

5.2 Altimeter And Attitude Control System Comparisons 

For the time period from 18 April 1975 to 30 June 1975, estimates of 
the GEOS-3 spacecraft attitude were • computed and published by the Goddard 
Space Flight Center [2]. These results were based upon magnetometer and 
sun sensor data from the GEOS-3 spacecraft and they provided estimates of 
the spacecraft z-axis pointing angle with respect to nadir. Due to the 
manner in which the altimeter antenna boresight was aligned with the space- 
craft z-axis, it is possible tha.t a 0.2 degree offset between the two axes 
exists [Private Communication, C. L. Purdy]. However, these spacecraft 
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ID HEADER: REV 184/022, COMPILE TAPE 7696 PILE 04, 06/02/76 


HHMMSS.SS 

ALT. ST AT. 

AV.PLAT. 

AV.ASG 

AV.RAGG 

DELTA 

PT.ANGLE 

AV.RTP 

SIGMA 

234832.15 

232.92 

.0129 

-.0172 

-72.808 

2.332 

-9999.99 

151.987 

-9999.99 

234842.39 

210.80 

.0563 

,0220 . 

-49.607 

.609 

-9999.99 

64.716 

-9999.99 

234852.63 

174.88 

.0605 

.0635 

-49.983 

-.050 

1.34 

62.890 

40.55 

234902.87 

160.16 

.0769 

.0645 

-51.197 

.161 

.83 

60.158 

39.02 

234913.11 

185.40 

.1072 

.1018 

-52.226 

.050 

1.12 

56.270 

43.45 

234923.35 

216.60 

.0028 

-.0211 

-83.004 

8.569 

-9999.99 

178.971 

-9999.99 

234933.59 

203.16 

.0482 

.0356 

-71.718 

.261 

.46 

64.840 

12.46 

234943,84 

204.00 

.0547 

.0444 

-72.822 

.189 

.75 

64.840 

12.32 


Figure 27. Sample listing from the ten second averaging program which computes the estimated 
pointing angle in degrees (VT, ANGLE) and cross section 0® (SIGMA) in dB. 



attitude estimates do provide data which can be used to approximately deter- 
mine the adequacy of the altimeter pointing angle estimation technique. In 
this section, a limited comparison of the two data sets will be made. 

5.2.1 Global Mode 

Figures 28-33 show comparisons between altimeter derived estimates and 
the attitude control system estimates for the pointing angle in Global Mode. 
Apart from the results in Figure 31, the altimeter estimates tend to always 
be larger than the spacecraft estimates. In particular, the altimeter de- 
rived estimate never goes below about 0.5°. While the data base is admittedly 
limited, the results in Figures 28-33 do tend to indicate that the altim- 
eter derived estimates are too large for a pointing error of less' than about 
0.6 in the Global Mode. It was initially thought that this discrepancy ^ 
might be due to a misalignment between the altimeter antenna boresight anH 
the spacecraft z-axis. However, subsequent inspection of Intensive Mode 
estimates of pointing angle did not corroborate this hypothesis. That is, 
the discrepancy appeared to be unique to the Global Mode of the altimeter. 

For the Global Mode Ramp and Plateau gates, gain calibrations were 
only obtained at room ambient temperature and pressure during preflight 
testing/calibration of the altimeter. Thus, the temperature dependent trans- 
fer gain of the Plateau gate xras unknown. Since this quantity has a direct 
Influence upon the estimation function, see equation (85), it was felt that 
inadequate calibration data could be a source of error. The analysis pre- 
sented in section 2.4 clearly shows that the estimation of pointing angle 
using (85) is relatively sensitive to changes in gain. In particular, the 
error is a nonlinear function of the pointing, angle; that is, for a 5% error 
in gain, the resulting pointing angle error is 0.4® at 0® pointing angle a nd 
0.1® at 1.0® pointing angle (see Figure 18). 

For the Global Mode, it is concluded that the pointing angle estima- 
tion technique using altimeter data is not accurate for pointing angles of 
less than 0.6® because of inadequate calibration of the Plateau gate gain, 
^ile it is possible to obtain a correction curve for the data (using In- 
tensive and Global mode comparisons) , this would entail a rather extensive 
effort and, hence, a reexamination of the need for Global Mode pointing 
angles estimates. 
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C (deg) 



Figure 28. Comparison of attitude control system derived pointing angles (A-^Ar-A) 
vith altimeter derived estimates (o— o— o) for Global Mode, Rev 167 
Start Time = 19H11M45S, 21 Apr 75. Standard deviation of altimeter 
estimate was 0.0506°. 
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Figure 29 . Coiaparison of attitude control system derived pointing 

angles (A-A— A) with altimeter derived estimates (o— o— o) 







. ■ Figure ’33. 
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Comparison of attitude control system derived pointing angle (A— &-A) with altimeter 
estimates ,(o-o-o) for Global Mode, Rev 416, Start Time = 09H06M25.84S, 9 May 75. 


5.2.2 Intensive Mode 


Orbit number 2762 was first selected for study due to the availability 
of sufficient data. Although the ten second averaging program had not as 
yet been developed, a GAP (version 2,0) Summary along with the GAP' and 
CAMIMERGE dximps were obtained. The CALIMERGE dump was used since it lists 
raw, uncalibrated data which can be hand-converted to account for the trans- 
fer gain of the integrating gates. The CALIMERGE results .were then compared 
with the GAP (version 2.0) Summary and GAP dump data. It was discovered 
that the GAP Summary incorrectly reported the Plateau gate value as 0.13 
input volts during the interval under consideration. Both the GAP dump and 
the hand converted CALIMERGE values agreed and indicated a lower value for 
per frame averages of .09 to .10 volts for the Plateau gate. It was appar- 
ent that there was a calibration problem in the GAP (version 2.0) Summary. 
This problem was further reflected in the failure of the pointing angle al- 
gorithm since for APG = 0.13 volts, the resulting value of A was larger than 
0.504 (see Table V). 

During this time, a revision of Reference [17] was obtained and while 
working with the I-Mode Plateau gate calibration table, an incorrect table 
entry was discovered. Under the 0®C temperature coltmm corresponding to a 
functional unit value of 0.1 volt, an engineering Tinit value of —0.99 volt 
is listed. This seemed unreasonable based upon the trend in values sur- 
rounding this entry and in fact appeared as though +0.99 volt ought to be 
the correct value instead. Using the CALIMERGE dump again for orbit 2762, 
a hand conversion from engineering to functional units was repeated using 
the revision for the Plateau gate calibration table. The results then 
agreed with the incorrect GAP (version 2.0) Summary values. Evidently, the 
GAP Summary had been using this incorrect table and consequently generated 
the xnrong Plateau gate per frame averages. Confirmation was made that this 
was an erroneous value and xmfortunately was being used in recent data 
processing. This meant that I-mode pass data having Plateau gate engineer- 
ing units and Intensive Tracker Temperature (ITT) values within range of the 
incorrect table entry were processed incorrectly for the conversion to Pla- 
teau gate functional units. 

Dr. G. S. Hayne of Applied Science Associates has developed a Fortran 
subroutine capable of taking Intensive Mode APG values, which have been 
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incorrectly calibrated for tracker temperatures under 20®C, and reversing 
the calibration procedure to retrieve the original engineering unit value. 
This value is then re-calibrated vith a procedure using the correct table 
value to output properly calibrated APG functional units. At present, this 
routine has been adapted to run with the ten second pointing angle program. 
Figure 34 is a typical program output before correction while Figure 35 is 
the same run after correction. It should be noted that the pointing angle 
correction also has an effect on the computed values for 

Despite the problem with the Plateau gate values, a number of passes 
were found to be operating in a temperature range such that the calibration 
was correct and the pointing angles derived from the altimeter data could 
be used to compare with attitude control system results. Figures 36 and 37 
illustrate typical comparisons of the altimeter and attitude control system 
results. Although the data is very limited, there is no apparent bias as 
in the case of the Global Mode. 

6.0 PROCESSING FOR AND INITIAL MEASUKEMENTS OF 

Closely associated with the computation of altimeter pointing angle is 
the estimation of the radar cross section per unit area, Cf®. Although the 
pointing angle is only computed every ten seconds, (T® can be computed more 
frequently. That is, once estimates of the pointing angle are available, 
it is possible to go back into the data and compute (J® at a much greater 
rate than one per ten seconds. However, the programs presently in use at 
WFC only compute O® at a ten second rate. 

6.1 Algorithm Development 

As shown by equations (62) and (64) , the average received power as 
seen by the Plateau gate can be written in the following form; 

For® (^ )1 

<Pj.(t)> = P^l^— j;-^Jf(?,Mode) (88) 

where the function f(^,Mode) depends upon both the pointing angle and the 
operating mode. A plot of f (?,Mode) was presented in Figure 19 for both 
the Global and Intensive modes. A closed form approximating function was 
fitted to the curves in Figure 19 in order to avoid the need for a look-up 
fable. Solving (88) for a®(i{;^)/L resulted in the following algorithms; for 
Global Mode. 
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ID HEADER; 

REV 1164/101, COMPILE TAPE 7169, FILE 

07, GETTAPE, 05/25/76 




HHMMSS.SS 

ALT.STAT. 

AV.PLAT. 

AV.AS6 

AV.RAGC 

DELTA 

PT .ANGLE 

AV.RTP 

SIGMA 


THE GATE 

VALUES WENT 

OUT OF RANGE 






63552.30 

79.00 

.0784 

.0509 

-70.075 

.351 

.71 

62.214 

11.84 

63602.54 

79.00 

.0756 

.0519 

-70.380 

.313 

.78 

62.214 

11,92 

63612.78 

79.00 

.0764 

.0524 ' 

-70.777 

.313 

.78 

62.214 

li.52 

63623.02 

79.00 

.0748 

.0525 

-70.781 

.299 

.81 

62.214 

11.67 

63633.26 

79.00 

.0746 

.0516 

-70.412 

.309 

.79 

62.214 

11.93 

63643.50 

79.00 

.0746 

.0515 

-70.114 

.310 

.79 

62.214 

12.22 

63653.74 

79.00 

,0732 

.0495 

-69.517 

.324 

.76 

62.214 

12.67 

“ 63703.98 

79.00 

.0731 

.0520 

-70.234 

.289 

.83 

62.214 

12.31 

63714.22 • 

79.00 

.0726 

.0513 

-69.601 

.293 

.82 

62.214 

12.90 

63724.46 

79.00 

.0729 

.0513 

-69.305 

.296 

.82 

62.214 

13.17 


Figure 34. 

Ten second 

pointing angle and o'* 

program output showing 

incorrect average 




Plateau Gate values (AV.PLAT.) re 

suiting from 

erroneous 

calibration table 

entry. 
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HHMMSS.SS 

ALT.STAT. 

AV.PLAT. 

AV.ASG 

AV.RAGC 

DELTA 

PT. ANGLE 

AV.RTP 

SIGMA 


THE GATE 

VALUES WENT 

OUT OF RANGE 






63552.30 

79.00 

.0952 

.0509 

-70.075 

.465 

.36 

62.214 

9.43 

63602.54 

79.00 

.0939 

.0519 

-70.380 

.447 

.44 

62.214 

9.34 

63612.78 

79.00 

.0953 

.0524 

-70.777 

.450 

.43 

62.214 

8.90 

63623.02 

79.00 

.0946 

.0525 

-70.781 

.445 

.44 

62.214 

8.95 

63633.26 

79.00 

.0949 

.0516 

-70.412 

.456 

.40 

62.214 

9.20 

63643.50 

79.00 

.0956 

.0515 

-70.114 

.461 

.38 

62.214 

9.44 

63653.74 

79.00 

.0951 

.0495 

-69.517 

.480 

.29 

62.214 

9.83 

63703.98 

79.00 

.0954 

.0520 

-70.234 

.455 

.41 

62.214 

9.39 

63714.22 

79.00 

.0953 

. .0513 

-69.601 

.462 

.38 

62.214 

9.95 

63724.46 

79.00 

.0961 

.0513 

-69.305 

.466 

.36 

62.214 

10.20 


Figure 35. Ten second pointing angle and 0 ° program output showing results 
of using the Plateau Gate calibration correction subroutine. 
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Figure 36. Comparison of attitude control system derived pointing angle (A-A-A) with altimeter 
estimates (o— o-o) for Intensive Mode, Rev 453, Start Time = 00H16M44S, 12 May 75. 
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Figure 37. Comparison of attitude control system derived pointing angle (A— A-A) with altimeter 
estimates (o— o— o) for Intensive Mode, Rev 530, Start Time = 00H27M00.4S, 17 May 75. 
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while for Intensive Mode, 


(89) 




(dB) = (AGC) 2 ^Qgg^ -10 log^g^l2.19 exp [- 2684.56 sin^^3 


cosh 


(3.201^ 


+ 152.15 - (RTF) 


lOsec 


(90) 


where the average in this case corresponds to 1000 pulses or ten seconds. 

In (89) and (90), C is in degrees, (^^^^lOsec transmit power in dBm 

averaged over a ten second interval, and the units of 0®/L are dB. The 
ten second average of the AGO, i.e., (^^^^Osec * be truly accom- 

plished in the following manner; the AGC output is converted from volts to 
dBm to remove AGC loop nonlinearities, the dBn values are converted to nu- 
meric, a ten second average of the numeric values is formed, and then this 
numeric average is converted to dan. In* the actual data processing, the 
dBm values of the AGC-inferred power are averaged since they do not exhibit 
a great deal of fluctuation. This same statement also holds for the trans- 
mit power, RTF. It should be noted, however, that for periods where the 
AGC output is rapidly fluctuating, the values of (O^/L^) resulting from (89) 
and (90) may be in error. For the Intensive Mode, the nominal angular lo- 
cation of the midpoint of the Flateau gate is ?{>^ = 0,28®, while for Global 
Mode it is i|^^-0.68“. 

For certain cases such as quick-look data reduction, it may be desir- 
able to decouple the 0® computation from the pointing angle estimation. 

That is, rather than averaging over a ten second interval to estimate the 
pointing angle, it may be more feasible to just list computed <7® values for 
some assumed pointing angle. This procedure woiild have the advantage of not 
requiring cross-frame averaging, which for quick-look data is certainly de- 
sirable. In this case, the assumption of 5 = 0“ might as well be made. Then, 
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for Global Mode 


gfXQ l ) . (dB) = AGC(dBiE) - RTF (dBm) + 149.14 (91) 

Lt 

p 

while for Intensive Mode, 

(dB) = AGC(dBm) - RTF (dBm) + 141.29 (92) 

Xj " 

P 

The AGC and transmit power outputs may either be averages or instantaneous 
values, depending upon the use of the data. Of course, this procedure is 
only intended for quick-look purposes where altimeter health or certain 
surface phenomenon require rapid examination of the data. Ignoring the 
pointing angle of the altimeter can result in more than a 3 dB error in cr“ 
for a one degree pointing error. 

6.2 Initial Cross Section Results 
6.2.1 Global Mode 

The results produced by the cross section algorithm used in the ten 
second averaging program were, for Global Mode operation, in close agree- 
ment with nominal values resulting from the Skylab data reduction effort 
[18]. Table VI gives one such comparison made for cross section calculations 
from SL-2 Pass 6 data on June 8, 1973, and the values obtained from GEOS-3, 
G-Mode Rev 415 on May -9, 1975, for Atlantic Ocean crossings. While these 
passes occurred two years apart, they are regionally similar and the sea 
states were not drastically different. These data may serve as a simple 
comparison of the values which, of course, should be system independent. 

It might not seem correct to present this comparison while an apparent bias 
in the Global Mode estimated pointing angle still exists. Justification is 
found, however, in the fact that the cross section computation is not very 
sensitive to small differences in the pointing angle up to about 1“ off- 
nadir. For instance, the first ten second average o'* value of 13.62 dB was 
computed for a 0.75" estimated pointing angle. With all other parameters 
held constant, a 0" pointing angle would result in a ten second cross section 
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TABLE VI 


Comparison of est im ated 0® values (in dB) taken 
from SKYLAB S-193 and GEOS-3 altimeter data, 
both from Atlantic Ocean crossings in Spring. 


SKYLAB SL-2, MODE V, 
PASS 6 
8 JUNE, 1973 

SM-0 

GEOS-3 
REV 
9 MAY 

G-MODE 

415 

1975 

FRAME 





NO. 

1-sec . 

HHMMSS.SS 

lO-sec. 

SSM-0 4 

12.7 

13.5 , 

073002.00 

13.62 

5 

12.8 

13.6 

073012.24 

13.85 

6 

13,2 

14.0 

073022.48 

14.69 . 




073032,72 

13.71 

SSM-1 1 

12.7 

13.5 

073042.96 

14.06 

2 

12.9 

13.7 

073053.20 

13.78 

3 

13,4 

14.2 

073103.44 

13.69 

4 

12.7 

13.5 

DATA BREAK — 

5 

12.5 

13.3 

073251.99 

13.24 




073302.23 

13.38 

SSM-2 1 

12.7 

13.5 

073312.47 

13.05 

2 

12.7 

13.5 

073322.71 

12.71 

3 

12.3 

13.1 

073332.95 

13.31 

4 

12.7 

13.5 

073343,19 

12.93 

5 

12.3 

13.1 

073353.43 

13.85 




073403.67 

13.61 




073413.91 

13.13 




073424.16 

12,45 




073434.40 

12.48 




073444.64 

11.92 




073454.88 

12.92 




073505.12 

14.36 




073515.36 

17.04 




— OUT OF 

RANGE — 

. 



073539.94 

13.68> 




073550.17 

12,66 


ROTE: The two columns for SKYLAB 0® values represent estimated 

bounds. The data rate for SKYLAB was 250 ptilses/sec 
and 1 Frame » 1 sec. 
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of 12.14 dB, less than 2 dB smaller. Given the degree of bias suspected In 
the estimated Global Mode pointing angle, one might expect the cross sec- 
tion computations to be about 1 dB too large. This 1 dB reduction places 
the GEOS-3 values even closer to the range of Skylab results in this case! 

6.2.2 Intensive Mode 

The initial application of the algorithm given by (90) to the Intensive 
Mode resulted in a® values which appeared to be more than 3 dB below Global 
Mode and Skylab results. For example, for relatively calm seas, the Inten- 
sive Mode values were below Skylab data acquired over 6 m. seas in the 
North Atlantic! Since the Intensive Mode estimates of pointing angle were 
in reasonable agreement with the attitude control system data, there was 
definite reason to suspect some other parameter in the algorithm. The use 
of the "Clean" AGC calibration curves* was suspected as the source of error 
and the theory presented in section 4 tended to substantiate this suspicion. 
When the "Clutter" AGC calibration curves were used to translate the AGC 
output into dBm, the resulting 0 ° values were in very good agreement with 
both Global Mode and Skylab data. The next section describes an in-orbit 
test that was designed to help resolve the "Clean" vs. "Clutter" AGC cali- 
bration problem. 

6.3 AGC Calibration Correction 
6.3.1 Rev 183-184 Crossing 

A fortunate situation for purposes of this study occurred on 22 April, 
1975, several hundred miles east of the coast of Iceland. Following Inten- 
sive Mode rev 183 (about 99 minutes later) Global Mode rev 184 crossed over 
rev 183’ s ground track and proceeded along for about twenty seconds in the 
same approximate area (see Figure 38) . This provided an interesting case 
for comparison of the two operating mode results based on the plausible as- 
sumption that surface conditions had not changed drastically in 99 minutes. 
In fact, a check on ground truth information [19] indicated moderate seas 


*See Section 4 for a discussion of the "Clean" vs. "Clutter" problem. 
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vTith 20 knot winds tinder cloudy skies and no weather extremes occurring in 
the vicinity of the crossing. It therefore seemed correct to expect cross 
section values derived from both passes to be reasonably -close in value . 

6.3.2 Cross Section Comparisons 

This section compares 0’“ values obtained from I-Mode (Rev 183) data 
with those obtained from G-Mode (Rev 184) in the vicinity of the crossing 
of each pass. Two different sets of I-Mode 0* values can be found using 
the two different AGC calibration curves in order to determine whether the 
".clean" test signal calibration or the "clutter" calibration will result in 
0® values which more nearly compare with those from G-Mode. As mentioned 
before, the close agreement between modes is expected by the assumption 
that cross section did not change considerably. In order to cope with the 
aforementioned problem in Global Mode pointing angle estimation, a pair of 
constant angle values were chosen to represent bounds for the actual values. 
For this purpose,, a lower bound of 0.2® and an upper bound of 0.7® were 
used. Then, 0® values were computed for each of the pointing angle values 
as a function of elapsed time along the pass. Table VII presents the re- 
sults of these calculations. Per frame 0® values were computed for each 
mode and in the case of the I-Mode calculations, the altimeter estimate of 
pointing angle for ten second intervals was used. As can be seen, the 0® 
results derived from Clutter curve converted AGC values are considerably 
closer to G-Mode values. Since the AGC conversion process using Clutter 
calibration curves results in I-Mode 0® values much closer to G-Mode re- 
sults, it is concluded that they should be used in data processing instead 
of the currently used Clean calibration curves. It would be well to further 
substantiate this conclusion with more comparisons using similar rev cross- 
ing situations. 

A "correction" routine has -been added to the ten second pointing an- 
gle and 0® program for use on the WFC ECLIPSE system. This routine accepts 
only Intensive Mode AGC values which have been processed using the Clean 
calibration data and produces a new AGC value which is equivalent to using 
Clutter calibration with the original raw AGC outputs. Because of the tem- 
perature dependence of this correction and the desire to keep it as simple 
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TABLE VII 


Comparison of Intensive and -Global Mode 0“ values during a segment from 
each pass crossing the same region of ocean. Intensive Mode values show 
the results of applying both Clean and Clutter AGC calibration cuarves. 


REV 183 REV 184 

INTENSIVE MODE GLOBAL MODE 


FRAME 

V 

'^clean 

clutter 

FRAME 

II 

• 

0 

o'»C5=.7“) 1 

77 

.33° 

5.51 

9.32 

89 

9.54 

10.76 

78 


6.12 

9.92 

90 

9.54 

10.76 

79 


5.86 

9.66 

91 

9.62 

10.84 





92 

9.54 

10.76 

80 

.24' 

6.25 

10.04 

93 

9.54 

10.76 

* 81 

- 

5.96 

9.75 

*94. 

9.94 

11.16 

82 


6.14 

9.93 

95 

10.14 

11.36 





96 

10.31 

11.53 

83 

.26° 

6.02 

9.82 

97 

10.18 

11.40 

84 


6.15 

9.95 

98 

10.31 

11.53 

85 


6.23 

10.03 

f 

99 

10.46 

11.68 




100 

10.65 

11.87 





101 

10.75 

11.97 


NOTE; The two data sets have been aligned as closely as possible. 
*These frames correspond to the point of intersection of Rev 183 and 
Rev 184. 
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as possible, the resultant values may, at times, differ from a straight- 
forward Clutter calibration by a dB or so. The additional time, effort, 
and computer storage required to reduce this error further still did not 
appear to be warranted. 
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APPENDIX 


The basic integrals to be evaluated are of the following form; 


OO 


(A-1) 


where n is an integer and u and ct are constants. To accomplish (A-l), it 
is first necessary to start with the case of n = 0, i.e. 


OO 

= J ' exp ^lix - a exp (x)] dx 


(A-2) 


and substitute 


Since 


aexp(x) » exp(z) 


X = - £a(a) + z 


and the limits are invariant under the transformation, (A-2) becomes 


I = exp 
o 


^-ViAn(a)] ^ exp Ijiz - exp (z)^ dz 


(A-3) 


The integral in (A-3) is given in [16] page 308 (equation 3.328); 

= exp [-vfn(a)] F(y) 

where F(*) is the gamma fimction. 

Combining equations (A-2) and (A-4) yields 
00 

J ' exp [iJx - a exp (x)] dx = exp [- y In (a)J F (y) 


(A-4) 


(A-5) 


Differentiating both sides of (A-5) n times with respect to y produces the 
following; 


1 nn 



CO 

expjjix -a exp(x)^ dx = -^-^lexp [-li £n(a)J r(ii)|^ (A-6) 

Equation (A-6) is the general result; for n-1 (A-6) becomes 

CO 

^ X exp [ux - a exp (x)] dx = exp £n(a)] Jlii(a)F(vi) + F' (y)| • (A-7) 


which is the result required in Section 4.2 for y=N and a“N/P^ . 
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